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Abstract 

A model of copper-oxygen bonding and anti-bonding bands with the most 
general two-body interactions allowable by symmetry is considered. The model 
has a continuous transition as a function of hole-density x and temperature 
T to a phase in which a current circulates in each unit cell. This phase pre- 
serves the translational symmetry of the lattice while breaking time-reversal 
invariance and the four-fold rotational symmetry. The product of time-reversal 
and four-fold rotation is preserved. The circulating current phase terminates 
at a critical point at x = Xc, T = 0. In the quantum-critical region about 
this point the logarithm of the frequency of the current fluctuations scales with 
their momentum. The microscopic basis for the marginal Fermi-liquid phene- 
menology and the observed long wavelength transport anomalies near x = Xc 
are derived from such fluctuations. The symmetry of the current fluctuations is 
such that the associated magnetic field fluctuations are absent at oxygen sites 
and have the correct form to explain the anomalous copper nuclear relaxation 
rate. Cross-overs to the Fermi-liquid phase on either side of Xc and the role of 
disorder are briefly considered. The current fluctuations promote superconduc- 
tive instability with a propensity towards "D-wave" symmetry or "extended 
S-wave" symmetry depending on details of the band-structure. Several experi- 
ments are proposed to test the theory. 
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1. Introduction 



Besides their exceptionally high superconducting transition temperatures, Copper- 
Oxide (Cu-0) based metals also have exceptional normal state properties.^ Landau 
Fermi-liquid theory^ and associated quasi particle concepts, which are a foundation 
stone for much of our understanding of phenomena in condensed matter appear to 
be inapplicable to their normal state. The principal problem is the development 
of a consistent theoretical framework in which the unusual metallic properties can 
be understood. Moreover, it is necessary to have a theory for the normal state to 
understand the mechanism for the superconductive instability. 

Every transport property in the normal state of Copper-Oxide metals has a tem- 
perature dependence unlike other metals. For example, the electrical resistivity has 
a linear temperature dependence down to for composition near the highest Tc for 
any class of Cu-0 compounds even when that Tc is as low as lOK.^ On the other 
hand, the equilibrium properties, such as specific heat Cy and magnetic susceptibility 
X are consistent with the usual temperature dependence and are in fact only about a 
factor of two enhanced over band structure calculations.^ The copper-oxide metals are 
thus qualitatively different from liquid ^He and heavy Fermion metals where strong 
interactions produce strong quantitative renormalizations in both equilibrium and 
transport properties without changing the asymptotic low temperature dependences 
and which are properly called Fermi-liquids. 

Recent re-examinations^ of the foundations of Landau's Fermi-liquid theory have 
confirmed the robustness of the theory for dimensions higher than 1 for any Hamil- 
tonian with non-singular low-energy interactions.^ Therefore a principal part of the 
theoretical task is to show that, in a model appropriate for Copper-Oxide metals. 
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elimination of the high energy degrees of freedom leads to a singular effective low 
energy Hamiltonian. The tremendous variety and number of experiments on Copper- 
Oxide metals severely constrain the form of such a low energy Hamiltonian. 

a. Constraints From Experiments 

A schematic generic phase diagram is drawn in Fig. (1) on the basis of the resistivity 
data. Where measurements are available every other transport property shows cor- 
responding regions. The insulating- antiferromagnetic phase near 1/2 filling and the 
superconducting phase are shown in bold lines. The normal state is roughly divided 
into four regions with dashed lines representing cross-overs from one characteristic 
temperature dependence in transport properties to another. Region 1 the non-Fermi- 
liquid phase has p(T) fa po + PiT and similarly remarkable "simple" anomalies in 
all the other transport properties. In Bi 2201 resistivity measurements^ are avail- 
able from lOK to 800K and in Lai.85Sr.i5Cu04 from 40K to 800K.^ The measured 
resistivity exponent in both is 1.05 ± 0.05. We may safely assume that it is 1 with 
possible logarithmic corrections. Such a behavior is observed only in a very narrow 
region near Xc. p(T) begins to decrease below the linear extrapolation as temperature 
is decreased in region 3 and is consistent with an asymptotic T^ dependence char- 
acteristic of a Fermi-liquid. On the low doping side,in region 2, there is a crossover 
to resistivity increasing with decreasing temperature. This regime may be termed 
insulating. It is almost certainly related to singular effects of impurity scattering in 
a non- Fermi-liquid (see Ref. (8) and Sec. (8) below). Strictly speaking, one should 
draw a third axis in Fig. (1) labelled disorder. If the zero-temperature intercept of 
the high temperature linear-resistivity is taken as a measure of disorder one concludes 
that it generally increases in the available data as x decreases. The limited systematic 

3 



data with independent variation of disorder and x clearly shows that impurities have 
a dramatic effect in the underdoped regime while in the overdoped regime, region 3, 
their effect is conventional.^ 

The size of the region 4 between region 1 and region 2 on the underdoped side 
depends upon disorder. The crossover between region 1 and region 4 is marked by 
a decrease in Cy/T and x with temperature with their ratio almost independent of 
temperature.^ In region 4 the resistivity drops below the linear extrapolation from 
region 1. This region has been termed the spin-gap region, but this is a misnomer. 
Not only do magnetic fluctuations, but optical conductivity and Raman scattering 
intensity in all measured polarizations decrease at low energies from their value in 
region 1. In region 2 this decline continues while the resisitivity begins to increase 
with decreasing temperature. 

^Prom data for YBa2Cu408(248) under pressure, ^'^ it appears that x is not a 
unique parameter for the crossovers in Fig. (1). The stochiometric compound 248 
has a Tc ~ 80K and a resistivity with a crossover from ~ T to a higher T dependence 
below about 200K, very similar to the properties of YBa2Cu306.7. Under pressure, Tc 
rises to llOK at P ~ lOOkbar. Simultaneously the resistivity becomes hnear down to 
Tc. The low energy excitations also change under pressure. In 248 the Raman spectra 
shows a low energy decrease. Under pressure as resistivity becomes linear the low 
energy spectrum is restored becoming like that of the optimally doped YBa2Cu306.9, 
i.e. frequency independent.^^ 

The schematic phase diagram of Fig. (1) suggests that the anomalous normal state 
(as well as superconductivity) are controlled by fluctuations around the point x f^i Xc 
and T 0. This is consistent with the marginal Fermi-liquid phenomenology^^ which 
suggests that the breakdown of Landau theory is due to scale invariant fluctuations 
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consistent with having a quantum critical point/^ QCP, (i.e. a singularity at T = 
0) near the ideal composition. The critical point itself is inaccessible due to the 
superconductive instability. The nature of the symmetry breaking at the critical 
point was not specified. 

The MFL has a single particle self-energy of the form, 



where x = max|a;| for |a;| >> T and = ttT for T >> |a;|, A is a coupling constant and 
cUc is a cut-off energy. The quasi-particle renormalization amplitude 



then vanishes logarithmically as {cu, T) — > 0. 

A microscopic theory should specify the nature of the critical point and the sym- 
metry on either side of it. It should also answer the question: If there is a critical 
point at X = Xc, T = 0, what about its continuation in the (x-T) plane? Shouldn't 
there be evidence of nonanalytic properties on a line in the (x-T) plane? Experimen- 
tally, there is indeed a cross-over in the properties in the (x-T) plane from Region 1 
to Region 4 of fig. (1). But why a cross-over rather than a transition? Or is it that 
the properties studied like transport and specific heat are often weakly sensitive to a 
transition? 

Equation (1) gives only the single particle scattering rate r^"^. This was used to 
understand^^'^^ the observed tunnehng conductance G(V) ~ |V| for T — > and to 
predict the lineshapes in single particle spectra. A crucial aspect of the properties 
of the Cu-0 metals is that the momentum transport scattering rates, measured in 
resistivity p(T) (as well as in optical conductivity^^ a{u!,T) and Raman cross-section^^ 



E(a;, q) = A ujin— + ixsgno; 



(1.1) 




(1.2) 
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SR(a;,T)) are also proportional to max (|cj|,T). So is the energy scattering rate 
measured by thermal conductivity «;(T).^® The experimental result that at (q, a;) 

r"^~r"^ r^T~'^r^T (1.3) 

sp mom en V / 

puts a strong constraint on theories. The single particle scattering rate is required 
in general to be at least as singular as the momentum scattering rate T^om- The 
immediate conclusion is that the experiments require Im X](<^, kp) ~ a;", a < 1. Angle 
resolved photoemission should be used to put stricter bounds on the single particle 
self-energy than have been done so far. But an easier way might be through the 
electronic spefic heat. The electronic specific heat is directly related to the exact 
single particle Green's function. For the marginal case, a = 1, i.e. Eqn. (1.1) gives 
Cv ~ N(0)T(1 + AWc/T). A more singular self energy, a < 1 gives Cv - N(0)T". In 
the experiments the electronic part of Cy is obtained only by subtracting the estimated 
phonon heat capacity and at x Ri Xc is reported to be consistent with ~ T. While 
logarithmic corrections to it cannot be ruled out, substantial singular departures are 
ruled out. 

The proportionality of the single particle and the transport scattering rates occur if 
the fluctuations leading to (1.1) are essentially momentum independent as suggested 
by the MFL phenemenology. (In that case there are no vertex corrections in the 
calculation of the conductivity. See further discussion in Sec. (7)). But this poses 
the serious dilemna that on the one hand we wish to be near a critical point, on the 
other that we need (nearly) momentum independent fluctuations. 

The trivial way to get (1.3) is if the experiments are in a temperature range 
T > a;* where uj* is the characteristic frequency of some fluctuations which scatter the 
Fermions. Then the density of such fluctuations is ~ T giving (1.3). This is ruled out 
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by the specific heat (and magnetic susceptibihty experiments) experiments.'^ If such 
fluctuations are of physical quantities like spin or density fluctuations of Fermions, 
a characteristic enhancement in 7 = Cy/T of 0(EF/a;*) ~ O(IO^) must occur as 
in the heavy Fermion compounds. Experimentally, the specific heat and magnetic 
susceptibility at the ideal composition arc consistent with Fcrmi-liquid behavior, ~ T 
and constant respectively, with no more than about a factor of 2 enhancements over 
non-interacting electrons in the measured temperature range. (The data however does 
allow for logarithmic or small power law corrections.) This as well as the fact that 
a{uj, T) and SR(a;, T) behaves smoothly in the range a; < T < a;c where cOc ~ 0(|eV) 
suggests that there is no low-energy scale near optimum doping and that the upper 
cut off frequency of the fluctuations is very high, 0(|eV). 

It is hard to imagine that the fluctuations due to the antiferromagnetic T = 
critical point at x .02 can have much to do^^ with phenomena at x Xc, with 
Xc > 0.15. Temperature independent magnetic correlation lengths of about 2 A are 
observed at x Ri Xc in YBa2Cu306.93^'^. Temperature-dependent lengths of ~ 20A at 
X fti Xc in Lai.85Sr.i5Cu04^^ but with less than ~ 10% of the total frequency integrated 
spectral weight in the q-dependent part. The normal state anomalies are identical in 
both compounds. Similarly, there is no evidence of any universal phase-separation 
fluctuations or charge density fluctuations^^ '^^ in different compounds at x fa Xc. It 
would appear that if a critical point is responsible for the unusual metallic state, it 
must be associated with some quite unusual order parameter which is hard to detect. 

Equations (1.1) and (1.3) cannot be used to understand the observed anomalies in 
NMR,^^ Hall effect, and magneto-resistance.^^ The anomalies in the Hall conductiv- 
ity, (Txy and magnetoresistance are quite revealing. In YBa2Cu306.9, where bandstruc- 
ture calculations^'' give a very small usual Hall conductivity due to particle-hole sym- 
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metry, the Hall angle ©h = o-xy/uxx, with magnetic field perpendicular to the plane, 
varies^^ approximately as between lOOK and 300K In the same range the normal- 
ized magnetoresistance Ap(H)/p varies^^ roughly as with O^/Ap/p f« 0(1). The 
Hall effect data^^^ in La2_xSrxCu04 however appears more comphcated where a sat- 
uration in the anomalous contribution occurs at temperatures below about 60 K. In 
view of this, it is not clear at the moment whether the Hall-effect anomaly is aleading 
low temperature singularity or an intermediate to high-temperature phenomena. It is 
worth noting that a singularity in the Hall angle ~ T~^, equivalently a Hall number 
diverging as T~^, imphes a spontaneous Hall effect in the hmit T — > 0, i.e. a Hall 
voltage in the absence of a magnetic field. Kotliar, et al.^^ have found this behavior 
in a solution of the Boltzmann equation which besides a scattering rate r^^^ ~ T, con- 
tains a phenomenological skew scattering rate proportional to the applied magnetic 
field which is ~ T~^. 

Perhaps the most astonishing of the normal state anomalies are the nuclear re- 
laxation rates ^urj-i ^^^^ ^Tjf^ of copper and oxygen nuclei respectively.^^ While 
(*""TiT)~^ appears to diverge as temperature is decreased suggesting singular local 
magnetic fluctuations at the O nuclei, (°TiT)~^ is a constant, which is the conven- 
tional behavior. Furthermore, (°TiTK)~^, where K is the measured Knight shift at 
oxygen, is within 20% a constant, irrespective of the compound studied or the density 
X in any given compound in the metallic range. Other experiments show that copper 
and oxygen orbitals are well hybridized. Nevertheless, it appears the local magnetic 
fluctuations at copper and oxygen sites are quite different. The conventional behavior 
on Oxygen, independent of x, rules out the scenario that the anomalies on Copper 
are due to the peaking of anti-ferromagnetic fluctuations. 

It is axiomatic that the fluctuations responsible for the anomalous metallic state 
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also are responsible for the instability to the superconducting state. The anoma- 
lous fluctuations develop a gap in the superconducting state as predicted^^'^'^^ and 
observed in a wide variety of experiments on quasiparticle relaxation rate deduced 
through transport experiments^^ and in angle resolved photoemission experiments.^^ 
The symmetry of the superconducting state appears to be consistent with ''D-wave" 
(if the lattice is assumed tetragonal). This issue is not completely settled yet.^^ 
Moreover the electron-doped material Nd2-xCexCu04 appears to be an "S-wave" 
superconductor . 

To summarize, existing experiments require an internally consistent microscopic 

theory to: 

(i) reproduce the phase-diagram of Fig. (1) with a non-fermi-liquid metallic phase 
near the composition for the highest Tc with cross-overs to fermi-liquid on the high 
doping side and with a strong tendency to insulating behavior due to disorder in 
the underdoped regime. The underdoped regime shows loss at low energies of both 
particle-hole excitations (in spin as well as charge channels) and of single particle 
excitations. 

(ii) have equilibrium properties like specific heat and magnetic susceptibility near the 
ideal composition consistent in the measured range of T with characteristic Fermi- 
liquid behavior to within small corrections and in magnitude be within factors of 0(2) 
of those for non-interacting electrons. 

(iii) have long wave-length transport relaxation rates used to interpret electrical con- 
ductivity and thermal conductivity at the ideal composition that satisfy Eq. (3). At 
X = Xc the fluctuations leading to the anomalous transport should have no scale other 
than a cutoff of the 0(1/2 eV). 

(iv) The fluctuations should have a symmetry such that they produce singular local 
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magnetic fluctuations at copper nuclei to give the observed anomalies in the copper 
nuclei relaxation rate, but no singular local magnetic fluctuations at the oxygen nu- 
clei, 

(v) The fluctuations should be capable of producing a pairing instability of D-wave 

symmetry. 

There are of course many other special properties discovered in a subject in which 
0(5x10'^) papers have been published. But I regard the requirements listed above 
as the most basic or a least the irreducible minimum. The theory developed in this 
paper attempts to meet these requirements and suggest a few crucial experiments. 

b. Choice of a Model 

The choice of a model with which to do microscopic theory should be influenced by 
the fact that copper-oxide metals are unique. None of the thousands of transition 
metal compounds studied share their properties. The point of view taken here and 
elsewhere^^ is that the unique properties of Cu-0 metals arise from their unique 
chemistry in which ionic interactions play a crucial dynamical role. This point has 
been extensively discussed^^ and will only be briefly repeated here. The divalent 
transition metal oxides at 1/2-fiUing can be put on the diagram^^ in Fig. (2) in which 
one of the axes is the normalized local repulsion energy U on the transition metal 



U is the energy to convert two transition metal ions with formal charge state 2+ to 
one with formal charge state l"*" and the other to S"*", while Ex is the energy to transfer 



U/W = [e(TM)^+ + E(TM)^+ - 2E(TM)2+] /W 



(1.4) 



where W is the bandwidth. The other axis is the ionic energy 



EJW = E(TM^+0^-) - E(TM2+02-) /W 



(1.5) 
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charge from the ground state configuration of a transition metal ion with charge 2"*" 
and a nearest neighbor oxygen ion with 2~ to transition metal ion with charge l"*" 
and oxygen ion with charge 1~ . Screening and dipole corrections, etc. in the solid 
are included in the definitions of U and Ex. 

As one moves from the left to the right of the periodic tabic, the ionization en- 
ergy of the TM falls thereby decreasing Ex/W relative to U/W with corresponding 
movement on Fig. (2). In the insulating state of Cu-0, Ex is only about leV. It is 
a charge-transfer insulator with the lowest energy one particle spectra primarily on 
copper, Cu^+ — > Cu+ while the one-hole spectra is primarily in oxygen, — > 0^~; 
see Figs. (3a and 3b) where the contrast to transition metal oxides towards the left of 
the periodic table is also shown. In the metallic state, obtained by doping, there are 
charge fiuctuations on copper and on oxygen of similar magnitude. Ex ~ leV is made 
up of from two sets of energies, the atomic level energies and the Madelung or ionic 
energies each of which is O(lOeV). Indeed all transition metal oxides owe their struc- 
tural stability to the ionic energies. But in the metallic state these energies have little 
dynamical role to play in most TM-oxides because there are hardly any fiuctuations 
on oxygen, i.e. although the ionic fiuctuation energy term in the Hamiltonian 

Ei<jVijninj (1.6) 

has J2i<j Vij ~ O(lOeV), the fluctuation 5noj on the oxygen ions require a large energy 
and are insignificant. An effective low energy Hamiltonian of the Hubbard form is 
then adequate. This is not true in the metallic state of Cu-0 where < Shq > / < 
no >~ 0(1). 

One of the aims of this paper is to show that finite-range interactions, if sufficiently 
strong, lead to qualitatively new features in the phase diagram of the model. The 
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one-dimensional version of the model has been extensively investigated by numerical 
methods.^^ Bosonization methods give incorrect results for the model for V's > 0(Ex) 
whereas as is well known the one-dimensional Hubbard model can be bosonized for 
any value of U. Therein lies a clue to understanding how in one and higher dimension 
the low-energy properties of the model are quite different for small V and for large 
V. Bosonization is inherently a weak-coupling method; it works in the (ID) Hubbard 
model because the properties at weak-coupling (small U/t) are similar to those at 
large-coupling. The numerical results show a change in properties as V is varied; 
for V > 0(Ex) they show the charge transfer instabihty described below, and growing 
superconducting correlation length as temperature is decreased, whereas Bosonization 
methods find V's to be irrelevant. (ID-models, do not have some essential features, 
discussed below, of the model in higher dimensions.) 

c. Preview of the Properties of the Model 

In Sec. (2) we discuss that it is enough to consider a two band model representing 
Cu and O bonding -b and anti-bonding -a bands (but with 0-0 hopping included) as 
illustrated in Fig. (4). In the hole representation, the chemical potential // is in the 
lower band as shown. For non-interacting electrons /i would be in the middle of band 
a at 1/2 fiUing as in La2Cu04 or YBa2Cu306. For hole doping, as for most Cu-0 
compounds, fi rises with doping x, for example in La2_xSrxCu04 or YBa2Cu306+x- 
We will find it important to consider the most general form of two-body interactions 
allowed by symmetry in the space of these two bands. When the strength of the 
interactions is on the same scale as the overall electronic bandwidth every term has 
a crucial role to play. 

If 11 were in the gap between bands a and b, the model is identical to the Excitonic 
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Insulator problem,^" which has been extensively discussed. We investigate here the 
model with /i in one of the bands. This changes the problem substantively. 

It is well known that interactions completely alter the one-electron picture at 1/2- 
fiUing; a gap develops around the chemical potential and the relative amount of Cu 
and O character of the occupied and un-occupicd states is drastically altered, fi stays 
in the band significantly away from 1/2- filling but the relative Cu-0 character of the 
occupied and unoccupied states is again expected to be quite different from the one- 
electron picture. If no change in the lattice symmetry occurs due to the interactions, 
this is formally describable by a q = instability of the one-electron band structure 
to a state in which 

~ ReEk-^x(k) (a+ bk,) 7^ (1.7) 

Here J^x(k) is a form factor expressing the relative Cu and O character of states in 
bands a and b. We will see that Tx is closely related to the relative average charge 
in Cu and O orbitals. The mean-field free-energy as a function of Tx is shown in Fig. 
(5a). With Tx 7^ 0, the orbitals must be rehybridized leading to new bands a and (3 
of the same general form as a and b, as shown in Fig. (4). This is expected to occur 
for strong interactions at very high temperatures for all x of interest just as it does 
at X = 0. For the most general interactions, such a transition is of first order, as for a 
free-energy of the form shown in Fig. (5a). Previous investigations'*^'^^ of the model 
had focussed on this charge transfer instability. 

We show that the model also has an interesting second-order transition in the 
Ising class to a state in which 

Ty ~ Im Ek. -^y (k) ^ 0. (1.8) 

The mean-field free-energy, for a fixed Tx 7^ 0, as a function of Ty for various values 
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of a parameter p which is a function of x and T is shown in Fig. (5b). This transition 
therefore occurs on a hue Tc(x) in the x — T plane. We will be specially interested in 
the properties in the vicinity of the point x = Xc(0) where Tc = which we identify 
as the quantum critical point. 

A finite Ty provides an additional relative phase to the wavefunctions at Cu and 
the two O sites in a unit cell. We will see that the ground state with a finite Ty 
corresponds to a four-fold pattern of circulating current within a unit cell with all 
cells staying equivalent. This is illustrated in Fig. (6). We may call it the circulating 
current (CC) phase. Translational symmetry is preserved but time reversal-invariance 
and four fold rotational symmetry are broken. The product of the time-reversal and 
the four-fold rotaion are preserved. 

In connection with the excitonic insulator problem and the Hubbard model at 
1/2- filling circulating current phases have also been discussed, which also break trans- 
lational symmetry, and go by the names of orbital antiferromagnets^^ or staggered 
flux phases. There are two important differences - (i) Our use of a Cu-0 model with 
more than two atom per unit cell allows a q = transition to a circulating current 
phase, so that lattice translation symmetries are preserved. There is no change in the 
symmetry of the band-structure, (ii) We discuss such a phase in the metallic state. 
This leads to a very special nature of the collective fluctuations near the QCP due 
to scattering of the fluctuations by low-energy particle-hole excitations at the Fermi- 
surface. The problem of determining the spectrum of the fluctuations is the same as 
the absorption in a degenerate semi-conductor with flnite mass electrons and holes 
and in the limit that the interactions are much larger than the Fermi-energy. We show 
that the logarithm of the frequency of the fluctuations scales with their momentum 
- i.e the QCP has dynamical critical exponent Zd = oo. The fluctuations are thus 
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independent of momentum to a logarithmic accuracy. This extreme quantum hmit 
is essential to understand the observed behavior in long-wavelength -low frequency 
transport properties summarized by Eq.(1.3) and its finite frequency counterparts. 

The model is stated more completely in Sec. (2). One can study the model 
as interactions are increased from zero or take a strong-coupling point of view and 
consider corrections about an infinite value of the interactions. Both approaches 
lead to the same "intermediate energy scale" Hamiltonian whose analysis begins in 
Sec. (4). Sec. (3) is devoted to the strong-coupling expansion. The low energy 
Hamiltonian is derived and analyzed in Sec. (5) and (6) where it is shown that a 
systematic and controlled analysis is possible. The physical properties of the pure 
model are investigated in Sec. (7), and in Sec. (8) a beginning is made to consider the 
effect of impurities in the properties of the model. Impurities are strongly relevant 
near a — oo transition. For arbitrarily small concentration, they convert the 
hue Tc(x) to a cross-over. We also show that for arbitrarily small concentration of 
impurities in a nonFermi-liquid, the Fermi-surface withers away at low temperatures. 
The density of states at the chemical potential is zero and the resistivity is infinity 
as T — > (unless superconductivity intervenes). 

Not including the effect of disorder the schematic deduced phase-diagram of the 
model is shown in fig. (7). In region I, the properties are determined by quantum 
fluctuations and arc that of a marginal Fermi-liquid with a cross-over to a Fermi- 
liquid regime in region III. In region II Ty 7^ in the pure limit. This phase should 
have fermi-liquid properties at low temperatures in the pure limit but with different 
parameters from that of region III. The transition between regions I and II turns 
into a crossover at arbitrarily small concentration of impurities. At low temperatures 
impurities are expected to lead to a further crossover in region III to an insulating 
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regime with zero density of states at the chemical potential Antiferromagnetism 
phase near 1/2-fiUing and superconductivity phase near Xc which are also properties 
of the model are not shown. 

In Sec. (9), a beginning is made to study the pairing instability due to the 
exchange of the circulating current fluctuations and a propensity towards a "D-wave" 
superconducting instability is indicated. 

In Sec. (10) the shortcomings of the theoretical calculations as well as the un- 
explained features of the experiments are highlighted. Experiments are suggested to 
test several features of the theoretical proposal, the most important is to observe the 
current pattern of fig. (6) by polarized neutron or x-ray scattering. 

2. Model for Copper-Oxides 

The basic building block of the Cu-0 compounds is the elongated CuOe octahedra in 
which the planar-short bond oxygens are shared at the corners to produce a layered, 
anisotropic three dimensional structure. The interlayer kinetic energy depends on the 
details of the structure. In the least anisotropic compounds the inter-layer bandwidth 
is O(10~^) of the intralayer bandwidth. Since the properties of such compounds in 
the temperature region of the normal state are the same as of those with anisotropy 
ratio O(10~^), a two-dimensional model is appropriate for the essential physics. The 
basis set for the minimum Hamiltonian is the dx2_y2 orbital on the Cu-ions and the 
Px and py orbitals on the 0-ions - see Fig. (8). In this basis set the Hamiltonian is 
written as 

H = Ho H- Hi -h H2. (2.1) 



16 



Ho is the kinetic energy 

Ho = E -1^ (ndi - npxi - Hpyi) - //(ndi + np^i + npyj) 

+ I]i,o-tpd d^i (px,i+|,CT — Px,i-a,(T + Py,i+|,CT — Py,i,-f,(T) ^2 2) 

+ tpp Px,i+a,(T (Py,i+f ~ Py,i-f 

+ h.c. 

Here (i) sums over the unit cells in the plane, refers to the Cu-d-orbitals and 
P(x,y),i±f ,0- to the oxygen px or py orbitals, which are neighbors to the Cu in cell i at 
a distance a/2. The relative signs in the kinetic energy take into account the phases 
of the orbitals, as shown in Fig. (8). 

Hi is the short-range part of the interactions 

Hi = J2i,a Udndicrndia + UpHpiaHpia- (2.3) 

H2 is the long-range part of the electron-electron interactions and the exchange 
interaction important only for nearest neighbors 

H2 = Ei,jVyninj + ^^..^VxSi-Sj, (2.4) 

where Ui = Udi or npx,i, Upy^i as appropriate. One must include the long-range Vy ~ 
|Ri — Rj|~^ for |Ri — Rjl 00 to keep the long wavelength charge oscillations at a 
finite plasma frequency. Otherwise, only the nearest neighbor Cu-0 and 0-0 need 
be considered. 

Diagonahzation of Hq gives the model "one-electron" band-structure, which has 
"Cu-0 bonding" and "anti-bonding bands" and a O — O "non-bonding band" . For 
tpp = 0, the bonding band -b and the antibonding band -a are decoupled from the 
non-bonding band, tpp 7^ is important for the results of this paper, but I consider 
only the bonding and the anti-bonding for simplicity and neglect the non-bonding 
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band. The eigenvalues and eigen-vectors of the bands kept must include the effects 
of tpp. It will be clear that smaller interactions are needed in the three-band model 
than in the simplified model for the important instabilities derived in Sec. 4. 

The dispersion of the bands a and b will be denoted by ea(k) and eb(k). These are 
sketched in Fig. (4). Expressions for them, derived pcrturbatively in tpp/tp^ arc given 
in the Appendix A. Their eigen-vectors are specified by the annihilation operators for 
states in them in terms of annihilation operators for d and px,y orbitals at lattice 
sites: 

aica = Uad(k)dko- + Uax(k) sin (^^^ P^^'^ Uay(k)sin j Pyko- (2.5a) 
bka = Ubd(k)dka + Ubx(k) sin j Pxka + Uby(k)sin j Pyka (2.5b) 

where 

Px(y)ka - ^ Z-i Pix(y)fT« 

Ri label the unit cells; the Cu atom in the unit cell i is taken to be at Ri. Expressions 
for these coefficients, also calculated pcrturbatively in tpp/tpd are given in Appendix 
A. 

Bandstructure calculations give tdp, 2tpp, Aq ~ O(leV); various spectroscopic 
methods give the local repulsion on Cu - U ~ 0(10 eV). As good an estimate of 
the nearest neighbor Cu-0 repulsion as any is 

Vnn>eV[e(leV)Ro], (2.6) 

where Rq is 1/2 the unit cell size and e(leV) is the measured long wavelength di- 
electric constant at an energy of 0(1 eV), which is fa 4. This gives the nearest 
neighbor Cu-0 interaction > 1.7eV. So, multiplied by the number of neighbors, 
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and considering polarization corrections etc., the characteristic Madelung energies 
controUing Cu-0 charge fluctuations is also of 0(10 eV). Calculations of Madelung 
energies^^ in the actual lattice support this estimate. 

The interaction energies and the overall electronic bandwidth are therefore of the 
same order. We can directly project the Hamiltonian (2.1) on to the basis set of 
the bonding and anti-bonding k-space orbitals obtained by diagonalizing (2.2). This 
results in the most general two-band Hamiltonian allowable by symmetry. Such a 
Hamiltonian (projected to interactions in the spin-singlet channel, which alone is 
important) is given in Eq. (4.1). The reader may at this point skip directly to Eq. 
(4.1) and the subsequent analysis. Although this approach is quite consistent for the 
metallic state, it is hard to derive the insulating phase near 1/2- filling from such a 
basis or to see that although the low energy physics at 1/ 2-fiUing of the general model 
is identical to the Hubbard model, it may not be so in the metaUic state. A basis of 
local real-space orbitals constructed in the strong-coupling limit to suppress some of 
the charge fluctuations is more convenient. The next section is devoted to deriving 
the projections of the Hamiltonian (2.1) to such a basis. The subsequent analysis, 
which is a simple generahzation of slave Boson methods^^ yields the Hamiltonian (4.1) 
as well. The results of this paper are qualitatively similar starting from either end in 
the ratio of interaction energy to the one-particle bandwidth. 

3. The Strong-Coupling Limit 

It is convenient to rewrite the Hamiltonian as an intracell part and an intercell part to 
calculate in the strong-coupling limit. We define a linear combination Dj^ of operators 
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on oxygen orbitals in a cell i which hybridize with the dx2_y2 orbital in the same cell: 

D£ = ;^(Pi+ax - Pi-ax + Pi+ay " Pi-aJ- (3-1) 

creates an orbital which also transforms as a dx2_y2 orbital about the center 
of the cell i. In the geometry of the Cu-0 lattice the orbitals created by Dj^ are 
not orthogonal for near-neighbor i. Wannier orbitals u!ia{r) can be defined which are 
orthogonal for different i. Such orbitals can be expressed as a linear combination of 
the orbitals created by Dj.. 
Ho can be re-expressed as 

Ho = Ho, cell + Hq^ 

inter— cell) 

(3.2) 

where because there is no hybridization at the zone center 

Ho, cell = J2i ^o(ndi - HDi). (3.3) 

We write the kinetic energy completely generally as 

Hq, inter-cell = ^if^Af^'ta^i'T + h.C.) (3.4) 

SO that it reproduces the one-electron bonding and anti-bonding bands. Here (/x, v) 
sum over d or D. As already discussed, use of non-orthogonal orbitals makes tij non- 
zero over a range larger than nearest neighbors. This conflict between the necessity of 
using non-orthogonal orbitals to handle strong local, but not just on-site, interactions 
and Bloch waves for a periodic lattice appears unavoidable. We need not dwell on 
this because the final projected Hamiltonian (4.1) for further analysis depends only 
on symmetry. 
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The interaction terms Hi, and H2, remain of the same general form in terms 
of Djo-'s as in Eq. (2.3) and Eq. (2.4) with redefinition of coefficients. We will 
simply regard that the interaction terms have been re-written in terms of D's, without 
changing the notation for the new coefficients. 

a. States in the Strong-coupling Limit 

The low-energy Hamiltonian for this model will now be derived. Since some of the 
important interactions are intracellular we specify a basis set in terms of the states 
of a cell by cutting off the kinetic-energy connection between cells and the long-range 
Coulomb interactions. 

Consider an average occupation of (1 +x) holes per unit cell as required by (1 + x) 
negative charges per unit cell assumed uniformly distributed by imposing a chemical 
potential /i. The minimum low energy basis must then include states with one-hole 
and with two-holes per unit cell. 

We define the zero-hole state (pQ^ |0 > as the closed-shell (spin zero) configuration 
in which the charge state of all oxygen ions is and of all copper ions in Cu"*". 
One-hole States: These are of two kinds: 

(i) di^^\0 >; a hole in the Cu-dx2_y2 orbital with energy Aq — /J,. Chemically, this 
is the spin 1/2 state Cu"*"*"© which is the nominal ground state configuration 
of the insulator. 

(ii) d2o-i|0 >: a hole in the orbital created by D[^, i.e. the oxygen-d-orbital, with 
energy — Aq — /i. (If such an orbital were localized on one-atom, the charge 
configuration would be Cu+0".) 
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Lowest Energy Two-hole States (f)t\0 >■ The lowest energy two-hole state is the spin- 
singlet state with one-hole in the Cu d-orbital and the other in the oxygen d-orbitals. 
The energy of this state is = V — 2//. (If the latter were confined to one-atom, 
the charge configuration would be Cu+0~.) 01*" should be thought of as a hard-core 

Boson operator. 

Neglected States: These are 

(i) Two-holes in the bonding combination of oxygen orbitals with energy — 2Ao — 
2/1 + Up. 

No essential physical difference arises if we include these two-hole states also in 
the low-energy sector. 

(ii) Triplet state with one-hole on Cu-d and the other on the oxygens. This is above 
0j^|O > by the exchange-energy, Vx which is O(Vij) for (ij) nearest neighbors. 

(iii) Two holes on the Cu-d orbitals with energy 2Ao — 2// -|- Ud. 

(iv) Three or higher number of holes per unit cell. 

(v) The zero hole state 0oi|O > with energy 0. 

At low energies, the allowed cells in a cell i must fulfill the completeness relation 
or constraint 

V^i+V'i + cj>t^i = 1 (3.5) 

where ipi = {d^ dii d2| d2j^)i. It is convenient to introduce Pauli matrices cr and r 
to specify respectively the spin (f J,) degree of freedom and the orbital (1 2) degree 
of freedom in the one-hole sector of the problem. 
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To derive a low-energy Hamiltonian, we must project (2.1) to states which fulfill 
the constraint (3.5). To this end, the bare operators di^. and Di^. are expressed in 
terms of the constrained operators through the identities 



div = 75 0i^d2i_^sgna (3.6) 
1 

V2 



X^t = ;)= 0[^dii_,sgn(7. (3.7) 



The intra- cell terms are transformed by noting that 

ndio- — + (3.8) 

nOi, = T^2ia + T^4>i/2, (3.9) 

where 

ni(T = d^"^ dia etc., and n^. = (3.10) 
b. Hamiltonian in the Strong-coupling Limit 

In terms of allowed states in the cell 

VondiUDi = Von^i, (3.11) 

and for i ^ j, 

VijUdiUDj = Vij(niin2j + n^in^j 
+niin<^j + n^in2j). 

The term J2i Vijniin2j is assumed already included in the definition of the difference 
Aq of the one hole states. Summing over i and j and using (3.5) the last two terms 
cancel the second term and renormalize E^. In general, the interactions of the di 
state and the d2 state with the neighbors are different and symmetry allows terms of 
the form 

II(ij)Vij(nii-n2i) n^j (3.13) 
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These renormalizc Aq downward proportionally to the density of two-hole states as 
found in Hartree-Fock and other previous calculations. ^^'^^ Equation (3.3) is simply, 

Ho,ceii = AoX^.(nn-n2i) (3.14) 

Consider now the intcr-ccll part of the kinetic energy. Starting from a configura- 
tion obeying the constraint, the inter-cell kinetic energy leads to configurations which 
preserve the constraint as well as those that do not. Consider first the former. These 
are necessarily processes which alter the two hole occupation in cell i and one hole 
occupation in cell j to one hole in i and two hole in j or vice- versa, for example 

<l>tdt^JO>^dtJ^\0>. (3.15) 

Therefore, Eq. (3.4) projected to the lower energy states gives 

Ho,i„tercell = IZ(ij),^,^^ ti^,j^(0i^d^^id+^j0j + h.C.) (3.16) 

The kinetic energy also operates on the one- hole states of neighboring cells i and j, 
creating disallowed states and the disallowed two- hole states on j. Eliminating 
such a kinetic energy term by a canonical transformation leads to an effective low 
energy interaction in the space of the allowed one hole states. This process is similar 
to that by which a Heisenberg exchange Hamiltonian is generated from the Hubbard 
Hamiltonian.^^ The new feature here is that the one-hole sector has a r degree of 
freedom as well as a <t degree of freedom. So there is an exchange in r space as well as 
cr space. The one-band Hubbard model produces an isotropic exchange Hamiltonian 
because the full Hamiltonian is invariant to spin rotations. This is again true in cr 
space here, but not in r space. The eigenvectors of the pseudo-spin r have in general 
different local energies, A 7^ in (3.14), and different transfer integrals. Also, there 
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are very many different intermediate states with no obvious rotational invariance. 
Here we derive the form of the effective interaction Hamiltonian from completely 
general considerations. An explicit derivation with calculation of the coefficients is 
given in the Appendix. 

The intracell kinetic energy terms in (3.4) which connect the allowed one hole 
states to the disallowed states may be written in terms of products of operators 

0Oidaai, 0^^jd«aj, a = 1,2 (3.17) 

where 0,^ create the zero-hole state and creates one of the disallowed two-hole 
states labelled by fi. The canonical transformation consists in eliminating the inter- 
mediate states 0oi02/ij (where now the allowed two-hole state 4>'l~ is included in /i). The 
most general pair-wise effective Hamiltonian is the sum over products of two kinetic 
energy operators with appropriate energy denominator. It has the general form 

Hint = IZa,r (^'^i^^'^j^^) ^a',T' (^^''^jVr'V'ia'r') (3.18) 

where {a =t, i) and r = (di,d2). In (2.23) J the sum over intermediate 

high energy (disallowed) states of matrix elements to such states divided by the 
corresponding energy denominators. 
Equation (3.18) can be rewritten as 

= - Eij Jij(l/4 - (T, ■ o-j)(A/4 - TiATj). (3.19) 

The Hamiltonian is isotropic in cr space and A expresses the anisotropy in r space: 

A^^, = x^x^/ (3.20) 
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In (3.19) 

T.i = (d+di-4d2)i, (3.21) 

Txi = (d+d2 + d+di)i (3.22) 

Tyi = i(d+d2 - d+di)i. (3.23) 

With axes defined as in (3.21)-(3.23) the most general form of A is such as to generate: 

Hint = Hxy + Hanis, (3.24) 

Hxy = [JzzTirj + J^(r;ri + h.c.)](l/4-o-i-o-j), (3.25) 
Hanis - (J.x(rM + ^x^i) + J±(T;Ti + TV:))(l/4-o-i-o-j). (3.26) 

The only conceivable terms missing in (3.24-3.26) are those linear in Ty. Note that 
Tyi = i (D^di — dj^Di). Therefore (remembering that Di and di refer to wavefunctions 
at different points in the unit cell i) Tyi represents a current distribution within the 
unit cell i. Terms linear in Tyi can not be generated from a time- reversal invariant 
Hamiltonian. Another way of seeing this is that if one has two bands as in fig. (4) 
whose states are created by linear combination of operators a;*" and b;'" respectively, 
the most general two-body interactions (with operators on sites i and j) are (ignoring 
spin) 

dbj dbj db[ , b j bj bj bi , bj bj 
bj cij bj B;j Bj bj bj 
a^aj^ajbi, b[*'bj'"bjai 

plus Hermitian conjugates of these. The terms in the first line can be rewritten in 
terms of t^t^, in the second line in terms of (t^tI -|- h.c.) and (t^t^_ -|- h.c.) and those 
in the third fine in terms of (T^iTxj + h.c.) just as in (3.24-3.26). 
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At 1/2-filling, x = 0, the state d^^i|0 > is not allowed in the low-energy subspace. 
Then ipi = (di-fdij^) only, and the low-energy Hamiltonian is obtained by dropping 
the T dependence in (2.24). The familiar Heisenberg Hamiltonian is then obtained. 
If one drops the r variable in the metallic state as well, the familiar t — J Hamilto- 
nian derivable from the Hubbard model in the strong-coupling limit is obtained. As 
discussed in the introduction, this is not justifiable for the parameters of the Cu-0 
problem. The r degrees of freedom make the problem richer and afford the possibility 
of new physics pursued in this paper. 

At this point it is useful to collect all the terms of the effective Hamiltonian 

H = Ei(Ao + Ej VijU^j) (nii - n2i) - Ei Ai(nii + n2i + n^i - I) - Ei(nii + n2i + n^i) 
+Ho jintercell + Ei E^n^i -I- Hint. 



Ai enforces the constraint and /i is introduced to fix the hole density at (1 -|- x). Hint 
is given by (3.24)-(3.26). 

c. Mean Field for the Slave Bosons 

We look for uniform mean-field solutions 



We also look for spin-singlet solutions in the bonds (i— j) favored by the kinetic energy 
term in (3.27) and the "RVB" decomposition.^^ 



(3.27) 



Ai -(Ai) 



A 



(3.28) 




(3.29) 



J (Ea,T Xr^'itr'^jar) (Ea'r' Xr'V'ja'T'V'ia'T') 
~ eij(EarXrV'iarV'jar) +eiy4Jij 



(3.30) 



where 



(3.31) 
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are mean-field amplitudes. 

The other decomposition of the interaction term in which mean-field amplitudes 
for Ti are introduced is more important to us. Note that given spin-singlets in the 
bonds (i — j), uniform spatial solutions in r-space are favored by the interactions in 
Eqs. (3.24), for Jj. > J^. 

We now diagonalize the bilinear terms in r space, i.e. the first five terms in Eq. 
(3.27) and transform to k-space. This introduces bands a and b. Let a^g., b^^ create 
particles in these bands: 

(3.32) 
(3.33) 
(3.34) 
(3.35) 
(3.36) 

Here tk is the lattice momentum transform of ty. The effective Hamiltonian projected 
to the bands a and b is given by Eq. (4.1). Wc denote the dispersion of the two bands 
due to this diagonalization also by eta and etb- They have the symmetry properties 
of the bandstructure in the one-electron approximation, i.e. Eqn. (2.5). 

4. Analysis of the Two Band Hamiltonian 

In the strong-coupling limit the intermediate energy scale Hamiltonian obtained from 
(3.27) after (3.28) and (3.29) is 

H = -(A + /i - 1| Eka{a'ta^ka + Ka^ka) - (A + /X - E^)n^ 
+ J2k,a ^kaa3']^a^ka + ^kha^ka^ka + Hint. 



^ka — 


Ukdika + Vkd2k(7 


bko- = 


-Vkdiko- + Ukd2k(T 


Vk/Uk = 


tan 


itan-^k/A 


tk = 


tk0^ 




and A = 


(A + V02) 
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Hint is given by 

Hint — Hxy + Hanis (4.2) 

(4.3) 



Hxy — X]k,k',q <^z(k, k', q)rzkqrzk'q 



+ Jj_ (k, k' , q) (r ^^qT ^^/q + TykqT yk'q) , 

^anis — Z]k,k',q <^zx(k, k', q)rxkqTzk'q + h.c. , , 

+ JX(k, k', q)(rxkqrxk'q - TykqTyk'q), 

where Tk,q are defined in a-b space, (not to be confused witli momentum transforms 
of (3.21)-(3.23)), 



''zkq — ^-i-q^-k t»k^qb 



k 



Txkq = i(4+qbk + bi^+qak) (4.5) 

''■ykq = "^(^k+qbk — bk+qak). 

As mentioned in Sec. 2, Eq. (4.1), follows directly from the "bare Hamiltonian" , 
Eq. (2.1) by transforming to the non-interacting bands using Eq. (2.5) (neglecting 
the terms in the first line, which do not play a crucial role) . The transformation from 
the bare U's and V's to J"'s using (2.5) is straightforward and not explicitly presented 
here. The only important point to note is that one should always include both the 
particle hole channels. Thus VudiUpj is written as 



V 
2 



^(d+di^)(p+,Pj<,0 + (d+pj^/)(p+,di^) 



+ one-electron terms (4.6) 



We note that the instabihties discussed here do not occur for the model with only on- 
site interactions, just as in the case of the Hamiltonian derived in the strong-coupling 
limit. In the strong-coupling limit the kinetic energy parameters t and A depend on 
the hole density x through the dependence of 0q on x derived below. This in turn 
makes the effective interactions J depend on x also. If (4.1) is considered directly 
derived from (2.1), the kinetic energy and the interactions are transformations of the 
bare terms. One can interpret the operators ako-, bko- as in (3.32) and (3.33) or as the 
bare bandstructure operators given by equation (2.5). 

29 



For short-range interactions, the jT's can be written as a sum over products of 
separable functions with the symmetry of the lattice. In terms of the leading such 
terms, we define. 



where rjX — (x, y, z). The form-factors J>,'s are obtained by Fourier transforming 
(3.26) and using the rotations (3.32-3.33) or directly from (2.1). 

Tjj's have a rather messy form. I assume that there is one particular lattice har- 
monic which dominates and henceforth drops the superscript s. Of course, the domi- 
nant harmonic can be determined only from a detailed calculation. Such calculations 
are not done in this paper, nor are they necessary for the principal calculations drawn. 

In (4.1) we have dropped the mean- field decomposition (3.31). It simply renor- 
malizes the kinetic energy in (4.1) without introducing any new qualitative features. 
The interactions term in (4.1) is purely in the spin-singlet channel. Due to lattice 
effects, coupling of the form TykqTzk'q and TykqTxk'q are also produced but they vanish 
as q — > and play no essential role. I have dropped such terms. 

a. Instabilities 

We look for instabilities in the bandstructure of the q = intracell excitonic nature 
due to Hint. They can arise only if the parameters in Hint are the same scale as the 
bandwidth. To calculate the properties of the new states, one introduces, as usual 
uniform (q = 0) mean field amplitudes, which will be determined variationally: 



J,cik, k', q) = ^ JrK^'.iK q)^^^(k', q). 



(4.7) 



s 



T = 




(4.8) 



2^xxEk«0 + ^k:o)-^x(k) 



(4.9) 
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Ty = -^:^yyEk«o-Tk:o)-^y(k) = EkTy(k), (4.10) 
where ^a(k) = ^a(k, 0). We also define the amphtude and angle 9 by 

T± = |T, + iTy|,tan^ = Ty/Tx. (4.11) 

The splitting of the bands always provides an effective field acting on r^. Therefore 
any interesting instabilities can only be in the — Ty plane. We therefore look for 
instabilities which determine the magnitude T^ of an order parameter in the Tx — Ty 
plane and its angle 9, with respect to the Tx-axis. 

If Hanis were ignored, the model is isotropic about the z-axis. There would be 
just one transition of second order nature with massless collective fluctuations. The 
coupling to Fermions of the collective modes would vanish in the long wavelength limit 
because these modes arise due to breaking a continuous symmetry. Hanis reduces the 
symmetry so that, as shown below, the general model has one first-order transition 
and two second-order transitions of the Ising variety as the parameters in the model, 
(including fi) arc varied. 

Before we proceed with the calculations, it is useful to discuss the excitation 
spectra for interaction strength less than necessary to cause the instability. Consider 
first only Hxy. The problem of the excitation spectra between a partially filled band 
(a) and an empty (or fully filled) band (b) has been investigated in degenerate semi- 
conductors^" and with the approximation of a dispcrsionlcss band (b) for the Fermi- 
edge singularities^°~^^ in the core spectra of metals. The absorption spectra is given 
by, see fig. (9a). 

X{cu, q) - E, Ek ^(k' q)Ga(k + q,u + u;)Gb(k, u) (4.12) 
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where A is the complete vertex in the particle-hole channel with energy-momentum 
(a;,q), and Ga and Gb are the (exact) single particle Green's functions. Using the 
fact that band b is empty (at T = 0), the sum over v can be explicitly carried out 
with the result 

X(^, q) ~ ^(^' q)Gb(k, v) |^=_a,+6ak-M ■ (4-13) 

For small interactions, there is a modification of the spectra at the threshold energy 
et = eb(kF) — A*- We are interested only in interactions large enough that an exci- 
tonic collective mode, which does not overlap the inter-band a-b transitions is pulled 
out. The simplest approximation for the calculations is to consider a rigid Fermi-sea 
which only serves to block out a part of the phase-space. This is the ladder diagram 
approximation for A, fig. (9b). In this case one obtains a sharp collective mode with 
spectral function ~ 5{ijj — a;ex(q))- This is a poor approximation for the lineshape.^^ 
The dressing of the exciton by low-energy particle-hole excitations at the Fermi- 
surface; the simplest processes are represented in fig. (9c); modifies the lineshape 
non-perturbatively. The problem has been solved exactly^^ in the recoil-less limit, i.e. 
for a dispersionless b-band where the interaction strength can be parametrized by a 
phase-shift 5(e). In this case A is also independent of momentum k and x is therefore 
independent of q. The absorption hneshape is as sketched in fig. (10a). Near the 
excitonic threshold it is given by 

X{uj) ^ {lo - idcxy^'"''^'^) ,u; > Ue^. (4.14) 

Here 5q is the phase-shift at the chemical potential, modulo tt - the phase-shift required 
to pull an exciton from the continuum. a;ex is determined by the details of bandstruc- 
ture, density of conduction electrons and the strength of the potential. Note that in 
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the weak-coupling limit the absorption lineshape has precisely the same exponent, 
but the absorption edge is at the energy et- Thus Eq. (4.14) may be regarded as 
the Fermi-surface singularity pulled down to cUex, or that the absorption displays the 
excitonic edge as modified by shake-off of low energy particle-hole excitations at the 
Fermi-surface. Therefore processes which smooth the Fermi-edge singularities will 
also smooth the excitonic edge. 

The effect of a finite hole mass or recoil^° on the Fermi-edge spectra is to smooth 
the singularity. Auger-processes now introduce a self-energy for Gb which is smooth 
on the scale of the recoil energy eb(kF) — eb (zone-boundary). All k's from the zone 
boundary to kp now contribute to the absorption for any given external q, momentum 
is conserved by particle-hole scattering on the Fermi-surf ace. Both Gb and A are now 
functions of k. The extra integrations in (4.12) then rounds off the singularity over the 
recoil energy. If interactions are strong enough to pull out an exciton, the excitonic 
edge must be similarly rounded off around the excitonic edge 6i;ex- This is shown in 
fig. (10b). Similar behavior must exist for a range of q from to order the difference 
from kp to the zone boundary. We may write qualitatively that 



where F is the smaller of the recoil energy or cuex- X (x) has the form (4.14) for x ^ 1 
but near x = 0, x(x) is a smooth function of x. 

The two important points in the above discussion are (i) that for large enough 
interactions, an excitonic state is pulled out with or without recoil, (ii) that recoil is 
always a relevant perturbation, smoothing out the singularity at the Fermi-edge and 
therefore at the excitonic edge if it exists. 

As the interaction strength in (4.1) increase a;ex(q) decreases. The bandstructure 




(4.15) 
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in Eq. (4.1) is unstable for interactions for which Re x(0, 0) oo. It appears 
difficult to get explicit closed form expressions for x(q, a;) taking into account the 
dressing of the exciton by low energy particle-hole pairs. In subsequent sections, I 
present explicit results with the frozen Fermi-sea approximation, and then discuss 
from general considerations the essential features of the exact x(a;, q). 

b. Anisotropies 

The interband interaction terms J^^^ and J'[^ play a quite different role than J'± and 
Jzz- It is convenient to first focus on these differences. The effective interaction 
as derived in the appendix is in general smaller than J'^x and J'yy. But the amplitude 
Tz is always finite because of the splitting between the bonding-b and anti-bonding-a 
band in (4.1) i.e. an effective uniform field ~ r^. The term proportional to in 
(4.4) is therefore approximated as 



Tz and Tx will therefore have no low energy dynamics. So, at low energies, the 
corrections to the mean- field approximation (4.16) are unimportant. The second and 
third term in (4.16) merely renormalize the bandstructure provided T^ and Tx are 
finite. The first term may be written as 



Consider next the term proportional to J'[^ in (4.4). In the mean- field approxima- 
tion we write it as 



[T,Tx + T, TxkJ^xz(k) + Tx T,kJ^xz(k) 



(4.16) 



Acos^, A^J,-^T,T^. 



(4.17) 



(:^!)-'[(Tx)^-(T,)^]. 



(4.18) 
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Therefore this term acts as a quadratic anisotropy field in the x — y plane if an 
excitonic state condenses. This anisotropy may be written as 

B(cos2^), B = (Ji)-^(T^)2. (4.19) 

In general, higher anisotropics are also generated from the starting Hamiltonian. 
I do not give their derivation but merely introduce, in the mean-field approximation 
a term, C cos^^ in the free-energy. The mean-field anisotropy free-energy is then 

Fanis = A cos ^ + B cos^^ + C cos^^. (4.20) 

A, B and C are functions of the density of holes x as well as, in general, of temperature. 

c. Condition for Instability 

(i) Rigid Fermi-sea As already mentioned only the ladder diagrams, fig. (8b) are 
considered in this approximation. 
The mean-field free-energy then is 

- (4.21) 

We will be interested especially in the vicinity of the hole density where the Fanis 
vanishes to leading order. It is convenient then to begin the analysis by ignoring Fanis- 

One thereby determines Tj^, T^ and (j). Fanis is then used to determine 9. 
The variational one-electron Green's function is 

Go(k,u;) = { ^-'^ \ = (a;-Ho(k))-i (4.22) 
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The mean-field band structure EkQ,^^ and the eigenvectors ato-, /^ko- are obtained as 
usual by diagonalizing Ho(k): 



(4.24) 



(4.24a) 



/ Qiko- \ / Ck Sk W ak<7 \ 

\ /3ka J \ -Sk Ck y V bka y 

Here 

Ck = cosAk/2, Sk = sin Ak/2, 
tanAk = J^±(k)T^/(eka-ekb + 2 J^.(k)). 

The leading term of ^±(k) = Sxy(k). 

Minimizing Fmf with respect to A and /i gives 

</'' + Ek,„.f(Ekm-/^) = l, (4.25) 

and 

(/)^ = x, (4.26) 

where f(z) is the Fermi-function. Equation (4.25-4.26) imply that the Luttinger the- 
orem on the volume enclosed by the Fermi-surface is satisfied. 

Minimizing Fmf with respect to T^, T_l and yields respectively 

^ + E.„f(E.™-.)^^0 (4.27) 
^ + i:..f(E.»-/.)i?f-0 (4.28) 
+ f + i:,„f(EK„-*.)^ = 0. (4.29) 

^From the appendix note that J'± ^ J^. We expect T^ to have only a minor effect 
which is determined mainly by the "external field" A. The stability of the mean-field 
approximation for 0i is ensured by a finite value for the "Boson" chemical potential 
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E<^ + Ao. In fact, apart from detailed quantitative issues, we need look only at equation 

(4.28) which can be re- written as 

^ + J^xd.jP ^'^-" - - g\ - = (4.30) 

The (approximate) condition for one electron band structure to be unstable is ob- 
tained by 

At T = 0, the right hand side of (4.31) is of order N(0) In (W + e^je^. Since the 
"threshold energy" et and the bandwidth W are of similar order, we need 2J71N(0) 
of 0(1) to have an instability. Wc will in fact assume that 2i7j_N(0) is large enough 
that the instability is at a very high temperature. Stability is achieved by Tj_ ^ 0, 
which corresponds simply to changing the relative Cu and O character of the occupied 
and unoccupied band. This result is no more than the statement that just as at 1/2 
filling the charge state away form 1/2 filling is determined by the electron-electron 
interactions, not just by the one-electron bandstructure. This behavior has been seen 
in a variety of earlier calculations. "^^'"^^ We will see below that when Fanis is considered, 
the transition to T^ ^ becomes of first order, except at two points. It is only near 
those two points that interesting properties can arise, 
(ii) Soft Fermi-sea 

Only the ladder diagrams are considered in the vertex A in deriving Eqs. (4.27)- 

(4.29) ; i.e. the Fermi-sea merely acts to restrict phase-space. As already discussed, 
this is a poor approximation. For sufficiently strong interactions the one-electron 
band structure is of course unstable even with the inclusion of low-energy particle- 
hole excitations, such as in Fig. (8c, d). But several details change. In general, the 
mean- field amplitudes Tx, Ty and Tz are functions of frequency. But as a variational 
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ansatz, Eqs. (4.8)-(4.10) may still be introduced. The considerations of anisotropy 
still continue to hold as in Sec (4.b). The mean-field free-energy (4.21) is an approxi- 
mation of the general case, where the T and (f) dependence (after integrating over the 
Fermions) may be written as 



Equations (4.27)-(4.29) are rigid Fermi-sea approximation of the general condition 



to determine the variational parameters T^, T^ and 0. 

The qualitative form of x(^? q) for interaction strength less than the critical value 
has been discussed in Sec. (4a) and illustrated in fig. (10b). Let us define p as 
the parameter (which is a function of the parameters in the Hamiltonian, (4.1)) such 
that the instabihty towards Tj^ 7^ occurs at a temperature Tc for p = pc(Tc). 
First consider T^ ^ 0. What docs the condition (4.33) for the instability imply for 
Im x{^i q) when the latter is overdamped and has the shape as in fig. (10b) rather 
than a delta-function as in the rigid Fermi-sea approximation. x(a;,q) has to satisfy 
the requirement 



For uj small compared to a;ex, Im x('^) 0) ~ uj while for uj large compared to a;ex, 
it is nearly a constant up to a cutoff uj^ on the scale of the Fermi-energy. Then by 
Kramers-Kronig transform the leading term in 




(4.32) 



det x~^{^ = 0, q = 0) = 



(4.33) 




(4.34) 
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As p ^ Pc(0), a;ex — >■ 0. Near this point, Im x(a;, 0) ~ sgn(u;). 

We see that recoil reduces the singularity as a; — > of x(a;, 0) near the transition 
at T = as p — > Pc(0) from the 5-function of the rigid sea approximation or the 
exact result (4.14) for the recoil less case. In terms of (4.14) recoil makes the phase- 
shift 6o at the chemical potential irrelevant. This appears quite unavoidable - on the 
one hand recoil cannot prevent the instability if interaction is large enough, on the 
other hand the singularity in (4.14) which is a Fermi-edge singularity is wiped out by 
recoil. The result is that the instabihty occurs with the least singular form possible: 
x(a;,0) ~ ln(a; -|- iO). As discussed in connection with Eq. (4.15), the smoothing of 
the excitonic edge occurs at any external q due to the mixing by Auger processes of 
inter-band particle-hole pairs over essentially the whole range of momenta. Thus for 
interaction energies large compared to the Fermi-energy as required for the instability 
and therefore also large compared to the recoil energy, the frequency dependence of 
the Im xi^y q) is nearly the same over the whole range of q. The imaginary part of 
x{uJ,q) is then reminiscent of the form of the Cooper pair- fluctuation propagator^^ 
above Tc which is ~ ia;/max(a;, T) over a range of q smaller than the coherence 
length i-G- the size of the Cooper pairs. Here, given the strong coupling required 
to engender the instability, the excitons have size of the order of the lattice spacing. 
So the ia;/max(ci;, T) form of damping is expected to persist over most of the Brillouin 
zone. 

The form for x(a;, q) near for T = as p ^ Pc(0) is thus 



x(^,q) 



1 CO , (jJ, 



-1 



+ In , , + + (Pc(0) - p) 



-1 



^max(a;,a;ex(p)) max(a;, a;ex(p)), 

^ (4.35) 

where ujc is an upper cut off energy, and k, provides the scale of dispersion. At a finite 
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temperature, we must use the fact that the u! and T dependence in x(a;, q) must scale 
as uj/T. So, for finite T > a; 

X{^,^, T) 



+ ln T^^^^A +«;V + (pe(T)-p) 



^max(T,a;ex(p)) max(T,a;ex(p)), 

(4.36) 
d. Determining 9 

We now consider the effects of Fanis- On minimizing Fanis with respect to 6, one finds 
that the equihbrium value ©m is given by 

Phase I : 0ni = for (A + 2B + 4C) < (4.37) 

and 

Phase I' : O^n = tt for (-A + 2B + 4C) < 0. (4.38) 

For C > 0, there occurs a second order transition of the Ising variety to 

Phase II : < < tt, for - A(2B + 4C < A. (4.39) 

© continuously rotates in phase II as A, B, C vary. 

Noting that A, B, C are in general functions of x and T, we may write the mean- 
field anisotropy energy as 

Fanis = Go(x,T)(^ - e^(x,t))2 + .... (4.40) 

where 

Go(x,T) = 0, for ±A + 2B + 4C = 0, i.e. at ei„ = 0,7r. (4.41) 

Therefore the transition at ©m = or tt is of second order. Let us denote the transition 
line (understanding that as x is varied we will be concerned only with either I to II 
or I' to II transition) by Xc(Tc). 
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Let us stay for definiteness in the vicinity of 9m = 0. For Gm = 0, the high 
temperature transition occurs as a first order transition with a real order parameter 
Tx- When ©m ^ 0, the mean-field order parameter is complex; Tx + i Ty. Ty 7^ 
implies that in the ground state, a current distribution occurs within each unit cell 
which has the same phase in every unit cell. 

The content of the mean-field theory is summarized in fig. (5), where the free- 
energy is shown as a function of Tx and Ty. Throughout the temperature region of 
interest Tx 7^ 0. As (p = ±A -|- 2B -|- 4C) is varied by varying x (and T), a second 
order Ising transition from Ty = to the circulating current phase Ty occurs. 
(We note in passing that we have arrived at a new class of statistical mechanical 
model for quantum-critical points.) 

We will discuss below how regimes 1 and 3 of the phase diagram of fig. (1) may 
be identified with the Phase I (or I') and regime 2 and 4 with the Phase II of the 
mean-field theory. The nature of the fiuctuations will be shown to vary in Phase I 
as a function of x and T leading to a cross-over in the properties from regime 1 to 
regime 3 in fig. (1). We will also show that the phase transitions between Phase I 
and Phase II also becomes a cross-over for arbitrarily small disorder. 

e. The Circulating Current Phase: 

©m 7^ (0, tt) implies that in the ground state, a current fiows in each cell. Since the 
momentum of the instability is zero, the current pattern respects lattice-translation 
symmetry. The current pattern within a cell can be deduced from the mean-field 
Hamiltonian with © 7^ 0, which is now Eqn. (4.23) with the substitution Tj^ — > 
Tj^e'®, Tj^e^'® in the off-diagonal terms. To calculate the current pattern, first find 
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the eigen- vectors of the conduction band with B 7^ 0: 

E„(k)dw = (cka + ^z(k)T,)aw + -Fx(k)Txe^®bw (4.42) 

Now use Eqs. (2.5) to express a, b in terms of dk and Px,y(k) using coefficients given 
in Appendix A. This yields 

ttka = Uad(k)dka + U„x(k)Pxk(7 + Uay(k)Pyk<7 (4.43) 

where u's are complex coefficients 

U„(d,x,y)(k) = R„(d,x,y)(k)e^^«^--)W. (4.44) 

There is no need to exhibit the comphcated expressions for these coefficients because 
the current pattern can be deduced from their general properties specified below. 

The current in a bond going from a copper site to an oxygen site in the x or y 
direction is 



jdx - f tpd Im (d,:^pi+x,a/ 

= f tpd Ek<kp cos Rctd RaxSin(0d(k) - 0x(k) ^ ' ' 



(The sum over k of the term with sin^ is zero using inversion symmetry.) Now note 
that tpd changes sign x — > — x and the sum in (4.45) is symmetric under inversion. 
Hence the current between a copper-orbital at a site i and an oxygen orbitals at 
i + |x and i — |x are equal and opposite. This holds also for the current between 
copper-orbitals at i and oxygen orbitals at i ± |y. 

The current between two oxygen orbitals in the same cell 

jxy = Il^(Pi+x,o-Pi+y,o-) .r.\ 

= ^ Ek<kp cos R6x(k) R«y(k) sin(0x(k) - 0y(k)) ^ ' ' 
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The sum in (4.46) is identical for the other three oxygen-oxygen bonds around a given 
copper atom, but as is evident from the phases shown in Fig. (8) tpp reverses sign 
cychcally in going around the four bonds, and therefore so does the current. 

The direction of the current between the copper and the oxygen orbitals and 
between the oxygen orbitals fixes the pattern show in Fig. (6). So together with 
breaking time-reversal symmetry, four-fold rotational symmetry is broken. But the 
product of the two is left invariant. 

Some further conclusions can be drawn from an examination of the u's. If tpp = 0, 
■0Q,x(k) = '0cjy(k). Then any (ipadO^) — '0ax,y(k)) can be removed by a unitary trans- 
formation without affecting the eigenvalues or the eigen-vectors. This is physically 
obvious from looking at Fig. (6); a current between copper and oxygen orbitals is 
meaningless in the absence of a current between the oxygen orbitals. 

We can also deduce that there is no contribution to the currents from states on 
the diagonals in the Brillouin-zone, ±kx — ±ky. Correspondingly, there is no change 
in the single-particle eigenvalues in the circulating current phase along the diagonals. 
On zone faces kx = or vr/a, ky = or 7r/a, the eigenvalues do change. The lowest 
lattice harmonic consistent with these symmetries is dx2_y2. So there is a change in 
the single particle spectra in the circulating current phase of dx2_y2 symmetry. 

The ground state current contribution of each state k depends on k and there is 
0(1) electron per unit cell in the conduction band. The orbital magnetic moment of 
the circulating current in each of the quadrants in Fig. (6) is O(0.05/xb), with the 
assumption that the average state contributes ~ 1/4//B per unit cell. 
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5. Collective Modes and Fermion-Boson Coupling 

We now consider the fluctuations in the state 7^ 0. They are interesting only near 
the I (or I') to II transition. There is always a finite effective field coupling linearly to 
Tz. So the fluctuations in the z-direction are always massive. The interesting modes 
are in the Tx — Ty space. So deflne 

5Tx,q = ^JxxEu«q + ^k:q)^x(k,q)-Tx (5.1) 

5Ty,q = -^J'yyEk«q-W-^y(k>q)-Ty. (5.2) 

The effective Hamiltonian determining the fluctuations is 

Hfluc = E.. {<. K) Ho(k) ( l'^^ ) + ^(<^T;^^5T,, + 5T;^^6T,,,) 

+ Ek,q->^±(kq)(a]^+qbk + b+^^ak) ((5T^,q + ST+_^) 
+iEk,q^±(kq)(a++qbk - b+^^ak) (5Ty,q + 5T+_^). 

+ Ha^is. (5.3) 

Again, let us ignore the effects of anisotropy to begin with but choose Gm = 0, 
i.e. T± = Tx. The spectrum of the fluctuations 

D°(q,a;) = {6TJT,){q,u;) (5.4) 
DO(q,u;) = {STySTy){q,u;) (5.5) 



is given in the frozen Fermi-sea approximation, (flg. 8b), by 

(5.6) 



DO-i(q,o;) \ =^ +Ek,jmq)P[Gaa(k + q,a; + i/)Gbb(k,i^) 



\ DO-i(q,u;) / ±Gab(k + q,a; + i/)Gba(k,i/)]. 

At a; = 0, q ^ 0, the equation for Dy is identical to Eqn. (4.30) determining Tx. So 
a long wavelength massless phase or current mode exists as is to be expected when 
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the anisotropy in the — Ty plane is zero. The poles of give the dispersion of 
the amplitude modes. Their frequencies near q are of order Tx; they will not be 
considered further. 

Let us now include the effect of anisotropy, but stay near I to II (or I' to II) in- 
stability. ^From (4.40), the anisotropy energy provides a quadratic term Go(x, T)(5Ty 
to the fluctuations. Including this effect 

^ uj^ — K^q^ + Go(x, T) 

where re/a is the order of Tx. The spectral weight of the collective mode m*/m is 



m 



m 



0(Tx/W). (5.8) 



Equations (5.6) are special cases, in the rigid Fermi-sea approximation of the 
general equation 

X-^(q,a;) = (5.9) 

which determines the fluctuation spectra, just as 4.27-4.29 are special cases of Eqn. 
(4.33) which determines the one-particle spectra through fixing Tx etc. In the frozen 
Fermi-sea approximation, there is no damping of the collective fluctuations - the ex- 
citonic resonances have a spectral function proportional to a 5-function. As discussed 
in Sec. (4a) for the case when Tx = 0, inclusion of low-energy particle-hole fluctu- 
ations change the spectral function of the excitonic collective mode in an essential 
way. The x(q, i^) including low-energy particle-hole fluctuations with Tx 7^ has the 
same functional form as discussed in (4c), but calculated with the new bandstructure, 
Eqns. (4.23), (4.24). The dispersion of the soft excitonic collective mode at q ^ 0, 
described by Dy(q, a;), again has the same form for G = as discussed in Sec. (4c), 
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leading to Equations (4.35) and (4.36). Including the effect of the anisotropy on the 
fluctuations, we have 

D^(q,c.) = Do^^ -c \r\..2„. ^ 



nTTP^ +^^ n ?T p J [ +«:V + Go(x,T) 

max(|a;|, T, Go) \ max(|a;|, T, Gq)/ I 



(5.10) 

Here Dq parameterizes the spectral weight of the fluctuation expected to be of O(Ep^). 
The second order transition occurs when Go(x, T) = as in Eqs. (4.37) - (4.39). 
(Henceforth Go is dimensionless, having been scaled by Dq^). The fluctuations have 
a finite frequency Gq above and a finite value below the transition (characteristic of 
transitions of the Ising class) . 

Equation (5.10) is crucial in the analysis below of the properties of the model. It 
is clear that the very singular result from the rigid Fermi-sea approximation, Eqn. 
(5.7) is quite incorrect. (It also gives properties in d = 2 which are too singular 
compared to experiment.) The combined effect of infrared processes at the Fermi- 
surface and recoil together with analjd;icity requirements has been discussed in Sec. 
(4c) to lead to (5.10). This justification is only heuristic. An evaluation of processes 
like in fig. (9c) to find A exactly appears very hard, if not impossible. Earlier, a 
three-body scattering approach to the problem was suggested. It might be possible 
to evaluate Dy(q, a;) systematically in such an approximation. Note that when 5q of 
Eq. (4.14) is zero, as argued here, Im Gb(c<j) ~ a;~^, at least for the recoilless case.^° 
This is consistent with the conjccture^"^ that a three-body resonance at the chemical 
potential may lead to the observed normal state anomalies. 

There are no other massless modes in the model. Earlier investigations^^'^^ of 
the charge transfer instability in the model found a diverging compressibility indi- 
cating phase-separation. If one adopts a short-range interaction model, the density 
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fluctuations have a dispersion u; ~ q. Near the critical point of the charge transfer 
instability, a low energy mode of (5Tx couples to such density fluctuations pushing 
the "electron sound velocity" to zero. In a model with Coulomb interactions, the 
density fluctuations are at the plasma frequency. Phase separation then does not 
occur (unless the system has inhomogencously distributed fixed (ionic) charges). 

Consider now the coupling of the Fcrmions to the low-energy collective modes in 
the vicinity of Xc(0) where a transition from = Tx to a complex order parameter 
Tx + iTy occurs. The coupling of the Fermions to the STy fluctuations comes from 
the fourth term in (5.3) and a similar term in (4.4). We must re-express the Fermion 
operators in terms of the low-energy Fermions created by a^^ by using the rotation 
(4.24). The coupling is written as 

Hf-b = Ek,q,. i g(k. q) <q,.«k,.(5Tyq + st;_^). (5.11) 

One finds from the fourth term in (5.3), (corresponding term from (4.4) introduces 
no important difference) 

g(k, q) = go •^.(k, q) sin |^^J^±S_Z^ j . (5.12) 

In (5.12) we have introduced go with dimensions of energy, so that STy is hence- 
forth dimensionless. go is expected to be of 0(Ef). In (5.12) J^y is given by (4.7). 
g(k, q) vanishes linearly with q. So even through we have an Ising transition, the 
coupling of the Fermion to the fluctuations vanishes at long wavelength. This occurs 
because the bands are split due to "external" fields and the fact that no TxTy or r^Ty 
coupling is allowed at long wavelengths. The vanishing of the coupling at long wave- 
lengths is however not fatal because in the fluctuation spectrum Eqn. (5.10), all q's 
are (to a logarithmic accuracy) equally important in properties which integrate the 
spectrum to energies of the order of temperature. 
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6. Analysis of the Low-Energy Hamiltonian 



The end result of the preceding two sections is a simple Hamiltonian for the calculation 

of low energy properties: 

H = Ek,a e(k)ai[^aka + D^-^q, u)ST+^STy^^ , . 

+ Ek,qig(k,qXq,«k,.((5Ty,q + (5T+_J. ^ " ' 

We now analyze the properties of this Hamiltonian in the regime near the I to H 

transition near T = 0. 

It is best to start by considering the simple physical processes depicted in figs, 
(lla-f). Here the wiggly lines denote the fluctuation propagation Dy(q, a;) and the 
solid lines the Fermion propagation Go(k, a;) — {aa'^) {k, lu) . 

The real part of 11(a) for q — > 0, u; — > renormahzes the mass G. This can 
be absorbed in a redefinition of Xp. Similarly the leading q dependence ~ q^ merely 
redefines the coefficient k in Eqn. (5.10). The imaginary part gives the usual Landau 
damping contribution ~ icj/vpq. If we use renormalized Fermion propagators in fig. 
(11a), the result is modified to 

ia;/max(Im T,{uj, T, G), vpq). 

We will show that in the non-Fermi-liquid regime Im E(a;,T) ~ max(a;,T). This is 
an additive correction to the imaginary part in (5.10) and is therefore unimportant 
at all q. This would not be true in the rigid Fermi-sea approximation in which the 
fluctuations in D°(q, a;) are completely undamped. (In that case D° has the same 
functional dependence on q and cu as the transverse electromagnetic field propagators 
in a metal.)^^ One can also examine higher order renormalizations, Fig. (lib) and 
fig. (11c) to conclude that they are irrelevant, the imaginary part of the process in 
(lib) is proportional to uj while (10c) is proportional to a;^. 
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Consider next the Fermion self-energy graph, fig. (lid). The imaginary part of 
the self-energy is easily seen to be 

ImE(q,cu) = N(0) /^d'^-^x r'deImD°(2kFX,cu-e)g2(2kFx) (6.2) 

Jo Jei 



e 

X 



, e u — e 

tanh — - + coth 



2T 2T 



where 



ei,2 = vpdql =F kpx) + Re S(ei,2). (6.3) 

It is found consistent to ignore E(e) in the right-hand side of (6.3). For T = 0, the 
limits of e integration are to a; (except for an ignorable region of x integration of 
0(a;/EF)). The x-integration, using (5.10) for Dq for Go = then leads to a constant, 
so that the final result is proportional to u sgn ui. Similarly for T ^ cj, the result 
is proportional to T sgn uj. The self-energy has negligible momentum dependence (if 
the Fermi-surface has no significant nesting) . These results are true for any dimension 
more than 1. The d-independence of the self-energy, and other properties in which 
these fluctuations are sampled over energies to the scale of the external frequency 
and temperature arises because in (5.10), the fluctuations are essentially local (even 
though Re D(0, 0) diverges as G~^). We may express both the real part and the 
imaginary part of the self-energy by an expression which interpolates between the 
a;/T < 1 and the a;/T > 1 limits: 



E(a;,T) = ttA 



(2a;/7r)ln(^^^^^^) +i7rT 



(6.4) 



Equation (6.4) may be useful in analyzing angle-resolved photoemission experiments 
discussed below. 

Consider next the vertex correction shown in graph (lOe). In the limit q — > flrst 



49 



and then a; 0, it is given by a Ward-identity (in the pure hmit), 

A'" = - (6.5) 

z 

where z is the quasiparticle renormahzation amphtude given by (1-2). In the "q-hmit" 
it is given by another Ward- identity in terms of dE/5k. Since E is found very weakly 
dependent on (k — kp), this is ignorable. 

For general cu and q, a finite vertex correction 0(g) non-singular as a function of 
uj and k — kp) is found. If the "bare" coupling constant g(k, q) is less than 0(1), 
this may be simply absorbed in the redefinition of g(k, q). One can formally devise 
^ schemes to keep such vertex corrections controlled. 

We briefly consider the renormahzation of Gq in Do(q, a;) due to anharmonic in- 
teractions. The leading contribution comes from 

u|5Tq,^n5Tq,,u;f (6.6) 

where u > is a phenomenological coefficient expected to be on the same scale as the 
upper cut off energy of the fluctuations. The self-energy of the modes proportional to 
u, fig. (llf) has the leading temperature correction proportion to uT (independent 
of d) for X Xc(0). This may be absorbed in Xc(T)) and suppresses the transition 
temperature. (Hence forth we can drop the superscript on G and D.) A proper 
analysis of the fluctuations near the transition line which changes from a quantum 
transition at T = to a classical Ising transition at high temperature has not been 
carried out. 

A very important general point to note in this connection is that the correlation 
length exponent at T = as a function of (x — Xc) is while it is 1 for the d = 2 
classical Ising model. This is expected to turn the transition hue to a cross-over for 
arbitrarily small disorder as discussed in Sec. (8). 
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We can estimate roughly the different regimes of fluctuations from the properties 
of the propagator D(q, a;). At this point no sophisticated analysis of the crossover 
between different regimes is attempted. Consider first the departure of the transition 
temperature of the circulating current phase from T = at x = Xc(0). We assume 
parameters are such that at T = 0, the CC phase occurs for x < Xc(0). This is 
provided by an estimate based on the parameters calculated in the appendix and the 
decrease of A with x. Then the transition temperature Tc(x), given by the divergence 
in Re D(0, 0) is at 

l^ln 

This gives 

Te ~ -G(T,) (6.8) 

i.e., the transition temperature is essentially proportional to — G(0). 

For finite Tc(G), in a narrow temperature ATc region near Tc(G), the fluctuations 
are characteristic of the classical Ising model. We have not investigated here how the 
width of this regime varies with Tc. In this regime, the classical thermal occupation 
of fluctuations (where |coth(c<j/2T) = (n(a;/T) + ^) ~ ^), determines the thermo- 
dynamic and other properties because the characteristic energy of the fluctuations is 
O(ATc) <^ T. If the characteristic energy of the fluctuations is much larger than T, 
the zero-point occupation of the fluctuations dominates the properties, (in this regime 
cotha;/T ^ 1). The physical properties in this regime are governed by the quantum 
fluctuations. Within this regime, we must distinguish when the characteristic scale 
of the fluctuations is given by temperature itself and when it is given by G(x). The 
former is the non-Fermi-liquid regime and the latter the Fermi-liquid regime. ^From 
the form of D(q, a;), the cross-over between the two occurs at T G(x) for G > 0. 



max(Tc, G) 



= -G for G < 0. 



(6.7) 
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The momentum integrated fluctuation spectrum for T ^ |G| gives a measure of 
the frequency distribution: 



Above the narrow critical regime near Tc(G), the properties are governed by the 
quantum fluctuations. Non-Fermi-hquid behavior is to be expected. For G < 0, i.e. 
the ordered side, the fluctuations have a gap. So Fermi-hquid behavior (but with 
unusuaUy smaU parameters) is to be expected in the pure hmit. We will discuss 
below that this regime is very sensitive to disorder. 

For G > 0, the fluctuations have a gap of 0(G) for T <C G. Fcrmi-liquid behavior 
is therefore to be expected but with parameters determined by G. 

The different regimes are shown in flg. (12). We are now ready to calculate the 
physical properties of (6.1). Given (6.1), all physical properties can be calculated in a 
controlled and systematic fashion because of the unimportance of vertex corrections 
in the Fermion-Boson scattering and the Boson propagator. For instance the sin- 
gle particle self-energy may be calculated self-consistently by using the renormalized 
Fermion propagator in (11a). The answer remains unaltered as in other problems 
with momentum independent self-energy. In the fluctuation spectra (5.10), all mo- 
menta arc equally important to logarithmic accuracy in the regime controlled by the 
quantum critical point: q scales as In a;. Formally this corresponds to a dynamical 
critical exponent za — > oo. This appears crucial to understand many of the observed 
anomalies in copper-oxide metals. It should be noted that the propagator (5.10) for 
G = is not the Zd = oo limit of the propagators discussed for example in Refs. (13), 
which are ~ {luj/c^ + q")~^ with z^ defined to be (a -|- P). 




(6.9) 



f In f for a; < T 
In ^i^, for a; > T. 
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7. Physical Properties 



The transport properties in regime I (see fig. 7) which are controlled by the quantum 
critical point and the cross-over to the customary behavior at low temperature for 
the overdoped case, region III, are calculated below. 

In the pure limit, region II should show Fermi-liquid properties but with different 
parameters from region III due to the alteration of states near the Fermi-surface by 
Ty 7^ 0. In the next section I argue that the transition between region I and II is 
only a cross-over and that at low temperatures region II is dominated by disorder 
such that the density of states at the chemical potential is zero. One should however 
expect a bump in Cy/T and x the I to II crossover. There is a decrease in low 
energy fluctuations in region II as Go in Eq. (5.10) is finite. But with the density of 
states at the chemical potential tending to zero at low temperatures due to disorder, 
a Fermi-liquid behavior may never be observable except in very pure samples. 

In the pure limit, although an order parameter develops in region II, it is by no 
means clear that there exist observable singularities in Cy/T or x the transition. 
Certainly at = 0, z = oo, no singularities exist. The crossover to Ising singularities 
(only logarithmic in the specific heat) at high temperatures for x far from Xc may 
occur with a very small amplitude at observable temperatures. This requires further 
work. 

a. Single Particle Spectra: Angle Resolved Photoemission Experiments 
(ARPES) and Single Particle Tunneling. 

A one-particle self-energy of the form (1.1) was suggested on phenomenological grounds 
While ARPES were soon found consistent^^ with this behavior, there has been since 
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then a considerable development in such experiments. Closer bounds should be put 
to this prediction. A useful formula to fit the self-energy which interpolates properly 
between a;/T <S 1 and uj/T ^ 1, while obeying analyticity requirements is given 
by Eq. (6.4). At low energies and low temperatures this behavior is modified in an 
interesting way by defects, as discussed later. 

Single particle tunneling has traditionally been a powerful tool for measuring the 
frequency dependence of the single particle spectra. Here, the interpretation of the 
tunneling spectra is complicated by the fact that if the self-energy is momentum 
independent its effect is not felt in the tunneling spectra unless the tunneling matrix 
element is momentum dependent. The situation has been amply discussed in Ref. 
(14) and need not be repeated. Under suitable conditions, the conductance as a 
function of voltage, G(V) — G(0) ~ Im S(V) ~ V as observed for tunneling in the 
c-direction. The observed G(V) varies weakly for tunnehng in the a-b plane. This has 
also been discussed. A new experimental development is the observation^^ by inelastic 
scanning tunneling microscopy that [G(V) — G(0)]/|V| increases as the distance of 
the tip to the surface increases thereby decreasing G(0). This is in accord with Ref. 
(14). 

As in the case of superconductivity through electron-phonon scattering, tunneling 
spectroscopy should serve to identify the spectra of the glue for superconductivity.^^ 
If the collective mode D(q, cu) is the glue, the tunneling conductance in the super- 
conducting state, G(V) — G(0) ~ |V| above the superconducting gap as observed 
in appropriate geometry. Quantitative verification of these ideas has been difficult 
because the slope of the conductance curve depends on "extraneous" factors as dis- 
cussed. But it should be possible in carefully designed experiments to normalize 
away the extraneous factors. After normalization the slope should depend only on 
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the coupling constant which determines the superconducting transition temperature 
T 

b. Long Wavelength Transport Properties 

As discussed in Sec. 1, an enormous constraint is put on a theory of copper-oxide 
metals by the fact that if the long wavelength transport properties are interpreted 
by kinetic theory or by semiclassical Boltzmann equations, the scattering rate for 
momentum loss measured in electrical resistivity and the scattering rate for energy loss 
measured in thermal conductivity have the same temperature dependence. Within 
experimental uncertainty, the single particle scattering rate measured in tunneling or 
ARPES also has the same temperature dependence. This is especially surprising in 
a theory in which the breakdown of Fermi-liquid theory is sought through a critical 
point where the long wavelength susceptibility diverges. One might imagine then that 
only long wavelength fluctuations or forward scatterings are important in scattering 
the Fermions, so that the backward scattering required in momentum transport would 
make the transport rate for momentum have a higher temperature dependence than 
in the energy transport. 

The backward scattering is enforced in momentum transport usually through con- 
sidering the two processes shown in (12). If for instance the Bosons are acoustic 
phonons, the leading T^ contributions to the d.c. resistivity of each of the processes 
in fig. (13) is exactly cancelled leading to a resistivity proportional to T^. 

The situation is quite different with D(q, a;) of the form (5.10). Then there is no 
cancellation to leading order between the self-energy and the vertex diagrams of fig. 
(12) because for energy transfer of order temperature, momentum transfer throughout 
the zone is important. If D(q, a;) were truly independent of momentum, as assumed 
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in the marginal Fermi-liquid phenomenology/^ (12b) would be identically zero due 
to the vector nature of the incoming and outgoing vertices. This is generally true for 
any "s-wave" scattering. With D(q, a;) of the form (5.10), the s-wave scattering part 
for a; T and k and k' both on the Fermi-surface 



is 0(1). One can show by an explicit calculation that the part of the process fig. 
(12b) for a; ^ T does not change the argument. The conductivity can therefore 
be calculated from figs. (12a) alone with just a numerical renormalization of the 
coefficient. Note that due to lattice effects, conservation of momentum with initial 
and final states on the Fermi-surface does not imply conservation of current. If 
for arbitrarily small cu, scattering occurs from a given state on the Fermi-surface 
to a substantial part of the Fermi-surface resistivity limited only by the density of 
fluctuation results. Since the density of states of the fluctuations is essentially a 
constant, a linear in T resistivity is to be expected. 

The calculation of electrical resistivity, optical conductivity, thermal conductivity 
and Raman scattering intensity from (6.1) is therefore essentially the same as done 
earlier, ^^^'^ with similar results apart from logarithmic corrections, 
(i) Optical Conductivity 

Optical conductivity as a function of u and T is a much more stringent test of the 
theory than p(T) alone. Such detailed comparisons with experiments have recently 
been done by Abrahams. Earher calculations were reported in Ref. (12b). For 
completeness and to show the quality of the flt to the experiments, the experimental 
results for the inplane conductivity deduced for untwinned single crystal is shown 
on the same scale with the calculations with indicated parameters in figs. (13a) and 




(7.1) 
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(13b), respectively. 

The microscopic theory from the strong-couphng Umit, Sec. (3) provides an ad- 
ditional important feature: The (intraband) optical conductivity sum rule is 

/ a{uj) duj = ujI (7.2) 

Given the constraint (3.5) the allowed density fluctTiations determining are only 
between the one-hole and the two-hole states 0i. ^From (4.26) the density of the 
two-hole states is x. Therefore u"^ ~ x(l — x). A proportionality of (7.2) to x for 
X < 0.2 has been noted experimentally.^^ 

In regime (111), where the integrated fluctuation spectra is ~ ^/G, a cross-over 
from p(T) ~ T to T^/G and a corresponding change in (t{uj) is predicted below T ~ G. 
In the presence of impurities D(q, a;) becomes optically active. The mid-infrared 
bump in the underdoped regime may be attributed to such processes. However, in 
compounds with chains, the chain conductivity^^ which appears to have a bump as 
well must be separated for meaningful comparison. 

(ii) Thermal Conductivity 

The graphs (13) with energy current external vertices give the thermal conductivity 
k(T). It follows the usual kinetic theory expression 

^ ^C.(T)(v|)r,h(T) (7.3) 

with T^{T) — AthT and Cv(T) a; T In ^. Ath departs from Amom by numerical 
factors due to the different angular averages in momentum and thermal transport. 
The Wiedemann Pranz ratio h{T)/T(t{T) is expected to be ~ In {uj^/T). 

(iii) Raman Scattering Intensity 

As has been discussed before the Raman intensity in a lattice has a part proportional 
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to the current-current correlation function and hence 



SR(a;,T) ~ (n(a;/T) + l)u;a{u;,T) 



(7.4) 



So a Raman intensity independent of frequency and temperature near the ideal com- 
position is expected as observed. Cross-over to a behavior linear in cu at low cu in 
regimes II and III is predicted and has been observed in regime (II) with a cross-over 
to regime (I) under pressure. In principle such SR(a;, T) is expected in all polariza- 
tions, the relative intensity may in general be quite different. 

The collective fluctuations (5.10) also couple directly in the Raman experiment. 
But since ImD(0,a;) ~ u!a{u!, T), this also gives the behavior of Eq. (7.4). 

c. NMR and Inelastic Magnetic Neutron Scattering 

The application of the theory to the NMR properties has been also described elsewhere, 
so only the principal point is summarized here. The current fluctuations (5.10) gen- 
erate an orbital magnetic field which vanishes at q = both at the copper site and 
the oxygen site as may be seen from fig. 10. At finite q, an orbital magnetic field 
proportional to q is generated at the copper sites but not at the oxygen sites. This 
is because around the four-fold co-ordinated Cu site a circulating current due to the 
electrons can be constructed to 0(q), but not at the two- fold co-ordinated oxygen 
site. This gives rise to an anomalous orbital contribution to the magnetic correlation 
functions at the Cu-sites 



where ao is the radius of the Cu d-orbital. The nuclear-relaxation rate calculated 
using (7.5) has the correct temperature dependence to fit the observations on Cu.^^ 
Oxygen nuclear relaxation rate follows the Koringa law. 




(7.5) 



58 



One of the most important aspect of the experimental results^"^'^^ is that the 
oxygen relaxation rate divided by the oxygen Knight shift does not vary either with 
X or from compound to compound within experimental uncertainty. Given this fact, 
and the fact that antiferromagnetic fluctuations, to the extent they are seen, change 
the position of their peak and their width with x, it is impossible to take seriously 
proposals which rely on the cancellation of such fluctuations at oxygen sites to account 
for the observations. A more robust symmetry is called for. In the picture presented 
here lattice symmetry never allows Xorb(q, at oxygen sites. 

The predicted XorblcL;*^); Eq. (7.5), can be measured by inelastic neutron scat- 
tering. Perhaps inelastic x-ray scattering can help distinguish the orbital magnetic 
fluctuations from spin fluctuations. Equation (7.5) predicts an unusually smooth 
q-dependence and scattering up to the high-energy cutoff at any q. The measured q- 
integrated magnetic fluctuation spectrum in Lai.85Sr.i5Cu04 is consistent^^ with (7.5). 
Further tests are suggested especially in compounds where nesting features of the 
bandstructure do not introduce sharp q-dependent features at low energies. A direct 
test of the theory would be inelastic neutron scattering experiments in several Bril- 
louin zones and transformation back to real space to deduce separately the magnetic 
fluctuations on oxygen and on Cu. Only a Fermi-liquid contribution ~ N(0)a;/qvF 
for cu < qvp and beyond should be seen on oxygen, which in appropriate range is 
negligible compared to (7.5), which should be seen only on Cu. 

8. Effects of Impurities 

The problem of disorder in a non- Fermi-liquid is complicated (and interesting). Only 
a preliminary treatment of some ideas is presented here to clarify aspects of the phase- 
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diagram of copper-oxides in the underdoped regime and the fate of the transition to 
the circulating current phase in the presence of disorder. 

First consider the effect of disorder in the transition from phase I to phase II, the 
circulating current phase. The Harris criteria^^ may be used in the classical regime 
of the transition to determine if quenched disorder, which varies the local transition 
temperature Tc(r) is relevant. This is derived by equating the free-energy contribution 
due to fluctuation in Tc(r) in a correlation volume to the pure fluctuation energy in 
the same volume. If 

di/ - 2 < (8.1) 

disorder is relevant. Here v is defined in terms of the correlation length as ^ ~ 
(T — Tc)~^ for a fixed x. In the Gaussian fluctuation regime, u — 1/2 while in the 
critical fluctuation regime for the d = 2 Ising model = 1, so disorder is relevant 

in the former and marginally irrelevant in the latter. This is expected to be true at 
asymptotically high temperatures far away from x = Xc(0) in the phase diagram, fig. 
(12). 

Now consider the transition at T = as a function of (x — Xc(0)). At T = only 
zero frequency fluctuations come to play. So the dynamical critical exponent cannot 
affect the relevance of disorder. The Harris criteria may be expected to therefore to 
be valid, but we should define u through (z/ = z/q), ^ ~ (x — Xc(0))~''°. Using (5.10) 
and noting that the argument of the logarithmic at a; = 0, T = log|x — Xc(0)| scales 
as q^, 1^0 — 0. The Harris criteria then suggests that disorder is strongly relevant. Not 
much definite appears to be known about the physical state when this is the case. 
The best guess is that the phase transition turns into a cross-over and that a glassy 
low temperature phase results with random local orientations of the order parameter. 
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This is quite reasonable when correlation lengths are short (z/q = 0); there is local 
ordering around each defect with no correlations building up between regions around 
different defects. 

To summarize the above, the correlation length exponent changes from to its 
classical Ising value at asymptotically high temperatures and large Xc(0) — x. Corre- 
spondingly one expects only a cross-over in the x — T plane to a glassy circulating 
current phase. 

The problem is even more interesting because the single particle excitations begin 
to acquire more singular self-energy than (1.1) or (6.4) due to defects. 

It was conjectured^ that a no n- Fermi liquid is an insulator for arbitrarily small 
disorder (resistivity — > oo as T — 0) if superconductivity does not intervene at a 
higher temperature. In recent calculations this conjecture has been supported by 
some systematic calculations.^^ The result of these calculations is that the impurity 
contribution to the resistivity in a marginal Fermi-liquid is proportional to inT below 
a cross-over temperature 



Here £o is the mean free path due to impurities calculated in the Born approximation. 
Below such an energy scale the density of states at the chemical potential also tends 

to zero, as {£nuj)^^. 

The observed low temperature resistivity in dirty samples, or in samples in which 
superconductivity is suppressed by a large magnetic field, (as well as the temperature 
dependence of the anisotropy in the resistivity) are consistent with these calculations.^^ 
These calculations, relay on using the marginal self-energy to calculate the impurity 
scattering vertex through a Ward-identity. The experiment gives a inT resistivity in 
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a wide range of doping in underdoped samples, where in the pure hmit the fluctua- 
tions have a gap and a marginal self-energy is not expected. One possible way this 
can happen is if the fluctuations acquire a finite low energy spectral weight due to 
disorder. This would be consistent with excitations in the glassy state conjectured 
above. 

9. Superconductive Instability 

It is only natural that the fluctuations responsible for the anomalous normal state 
also lead to the instability to superconductivity. We again look to the low-energy 
Hamiltonian, Eq. (6.1), to deduce the effective interaction in the particle-particle 
channel. As usual, this gives for total momentum of the pair equal to zero: 



Equation (9.1) is now used to deduce the symmetry channel with the largest pair- 
ing interaction. The procedure followed is the generalization to more than one atom 
per unit cell case of that in Ref. (65), where it was shown that antiferromagnetic 
fluctuations promote even-parity spin-singlet pairing which has "d-wave" symmetry 
in metals with appropriate band structures. The situation here is much more com- 
plicated; a preliminary analysis is given below. The propagator D(q, a;) is to a very 
good approximation independent of momentum for frequencies cu of importance for 
pairing which are always higher than T. So D can be regarded as a constant, Dq with 
an upper frequency cut off lJc- The effective pairing Kernel is then 



H, 



pair 



Ek,k'g(k,k') g*(-k,-k')D(k-k',a;) 



(9.1) 




g^Do|J-y(k)+ J-y(k'|W 




) 
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where the upper case corresponds to spin-singlet (even parity) and lower to spin-triplet 
(odd parity) pairing. We wish to express this as a sum over products of functions of 
k and k' which have the symmetry of the lattice and which are mutually orthogonal. 
Sticking to the lowest lattice harmonics, we look for coefficients in the expansion 

(3) [Js + JAA(k)A(k')+JDD(k)D(k') + ....], (9.3) 

(-1) [Jt(Tk,Tk,, + TkyTk'y) + ....] , (9.4) 

where Js is the coefficient for the simple s-wave pairing, Ja for pairing of "extended 
s-wave form" , i.e. 

A(k) = coskxa -|- coskya (9.5) 
Jd for pairing of "d-wave form" , i.e. 

D(k) = coskxa — coskya (9.6) 

and Jt for the odd parity form sinkxa or sinkya. The in (9.3-9.4) refer to higher 

lattice harmonics, i.e. periodic functions of 2 kxa, 2kya and so on which we ignore and 
which arc automatically mixed in to the gap function below Tc due to the nonlincarity 
in the gap equation, sin^ ^ ^k-A,,/ j 2^ jg ^^^.j comphcated function, but it has 
two properties which help write down the leading dependence on k, k' consistent with 
lattice symmetry. It is zero for k = k' and peaks when the difference momenta is 
maximum possible, i.e. at (kx = = ^, ky — = n/a.. The lowest lattice harmonic 
satisfying these conditions is 

(9.7) 



1 - - 
9 



cos(kx — k^)a + cos(ky — k' 
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In l^y(k) + ^y(k')p, the only part I keep is a constant; the others give harmonics. 

One then gets the relative magnitudes in units of ggDo ~ AN~^(0) 

s-wave pairing: 3 Js = 3 

D-wave pairing: 3Jd = — 3/2 

extend s-wave pairing: 3Ja = — 3/2 

triplet pairing: —1 Jt = +1. 

This immediately implies that simple s-wave pairing and triplet pairing Kernels are 

repulsive. In the present case simple s-wave pairing is disallowed simply from the face 

that the effective interaction vanishes at long wave length, the triplet is disallowed 

because the fluctuations conserve spin. The Kernel for D-wave and extended s-wave 

arc attractive and of equal magnitude. The situation is thus identical to the case of 

antiferromagnetic fluctuations with fluctuations peaking at (kxa = tt, kya = vr). 

The Tc is determined, as usual from the linear gap equation projected to the lattice 
harmonics. T^ for D-wave is in general different than T^, depending purely on the 
bandstructure, and the chemical potential, exactly as for the case of antiferromagnetic 
fluctuations. For that case and with Cu-0 bandstructure on a square lattice d-wave 
is formed to be favored in explicit calculations. The large density of states due to the 
proximity to van Hove singularities in the (tt, tt) direction favors d-wave. The same 
is therefore expected in the present case. Just as for antiferromagnetic fluctuations, 
variations in the bandstructure near the Fermi-surface gives superconducting states 
of different symmetries for the same interaction vertex. 

The upper cut-off ujc of D(q, a;) is of 0(Ef); from flt to normal state transport 
experiments, the coupling A ~ 0.5. As for normal state transport, vertex corrections 
are unimportant for calculations of Tp etc. So a consistent theory can be built. It 
should be noted that when D(q, cj) acquires a gap or the material has a propensity 
to insulating behavior, Tc must go down. Tc is therefore maximum near x = Xc(0). 

64 



It is worth noting that a signature of the glue for superconductivity is provided 
by the tunnehng conductance. Under appropriate experimental conditions, as dis- 
cussed in Sec. (7a) and Ref. (14), dG(V)/dV (above the superconducting gap) is 
proportional to the density of state of the glue for superconductivity weighted by the 
q-dcpcndence of the coupling constant: the famous "a^(cj)F(c(j)".^'^ The present the- 
ory predicts this to be a constant (with small corrections) up to the cut-off Uc- This 
is indeed observed but only in some geometries for reasons discussed in Ref. (14). 
Further systematic studies are called for. Optical conductivity in the superconducting 
phase for frequencies larger than twice the gap also can be used to deduce the glue 
for superconductivity.^^ The existing data is again consistent with the form (5.10). 

The results of Ref. (65) show that any Bosonic fluctuations such that the pairing 
interaction is miminal for q = and maximum for q = (tt, tt) produces d-wave pairing. 
To distinguish between mechanisms requires the experiments discussed above. 

10. Concluding Remarks, Further Theory, Further 
Experiments 

This investigation has been based on two basic assumptions: (i) Breakdown of Landau 
theory in more than one dimension requires scale-invariant low energy fluctuations, 
(ii) The solid state chemistry of Copper-Oxide is special and responsible for its special 
physical properties. Accordingly, I have formulated the Copper-Oxide model of Sec. 2 
and tried to investigate its properties in a systematic manner to find unusual singular 
low energy fluctuations. The model does have an antiferromagnetic instability at small 
doping. It probably has other flnite q instabilities at larger x for a range of parameters, 
especially if there is a nesting of the Fermi-surface. Given the experimental data, I 
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do not regard the singular fluctuations near such instabilities as a solution to the 
fundamental problems stated in Sec. 1. I have found a q = transition to an unusual 
Circulating Current Phase on a hue the x — T plane in the general model in the 
pure limit terminating at a quantum critical point at x = x^., T = 0. The model 
has unusual low energy fluctuations in which the logarithm of the frequency scales 
with the momentum, so that the fluctuations are essentially local in space. Such local 
fluctuations are essential to understand the peculiar normal state transport anomalies 
in which the momentum scattering rate, energy scattering rate and the single particle 
scattering rate are all proportional to T. They have the right energy scale to give a 
"continuous behavior" in optical and Raman conductivities from zero frequency to 
energies of O (1 eV). The current fluctuations also produce local orbital magnetic 
field fluctuations which have the symmetry and temperature dependence to account 
for the extraordinary NMR relaxation rates on Cu and O nuclei. The fluctuations 
also couple to Fermions to give a superconducting instability with d-wave symmetry 
favored. Disorder appears to convert the transition line to the circulating current 
phase to a cross-over line due to the quasi-local nature of the fluctuations. 

There exist several incomplete aspects of the theory presented here. While it has 
been shown conclusively that the circulating current instability does indeed occur in 
the model, the phase diagram in the T — x plane has not been determined. This 
requires an explicit numerical solution of the mean-field equations (4.27)-(4.29) with 
an assumed set of reasonable parameters. I have relied on Refs. (50-52) and general 
analyticity condition for fluctuations near an instability to present a heuristic dervia- 
tion of the form of the fluctuation spectrum, Eqs. (4.35)-(4.36). A better calculation 
is desirable. A detailed treatment of the different regimes of fluctuations and the ef- 
fect of disorder is needed. A complete examination of the superconductive instability 
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Tc(x) is possible and should be done. 

What are the principal experiments on which this paper has been silent? First 
is the question of the very interesting magneto-transport anomalies. ^^'^^ I have in- 
dicated in section 1 that the experimental results do not appear to show that they 
have asymptotic low cncrgy-tcmpcraturc singularities. But even so, the peculiar sub- 
leading behavior ought to be calculated. It is true that circulating current fluctuations 
lead to chiral scattering in a perpendicular magnetic field. As noted in Ref. (28), the 
temperature dependence of such chiral scattering is reflected in the magneto-transport 
anomalies. But so far I have not succeeded in formulating their effect consistently. 
Second, very interesting changes in angle resolved photoemission spectra®^ have been 
observed in going from region 1 to region 4 of the phase diagram of Fig.(l). It would 
be very natural to try to associate these with the transition or crossover to the cir- 
culating current phase. The observed changes in the spectra are most pronounced 
where the Fermi-surface of the ideally-doped samples crosses the (tt, 0) — (tt, tt) direc- 
tion and least pronounced where the Fermi-surface crosses the (0, 0) — (vr, 0) direction, 
i.e. the changes have x^ — y^ symmetry. As discussed in Sec.(4e) the changes in the 
one-particle spectra in the circulating current phase also do have x^ — y^ symmetry. 
This is quite intriguing but a calculation of the one-particle spectral function in the 
circulating current phase is necessary to draw any conclusions. Third, there are as- 
pects of NMR experiments, especially the anisotropy in the relaxation rate which are 
not explained in Ref. (61). Understanding anisotropy effects in NMR requires a theory 
of the coupling of fluctuation between different planes. 

Are there experiments left to do after the 0(5 x 10^) already published to test 
the conclusions of this paper? The answer is yes, but most of them are difficult 
experiments. 

67 



The most direct and convincing test of the theory would be the observation of 
the circulating current phase in the circulating current phase in the under-doped 
samples and its evolution as a function of temperature. As discussed earlier, long- 
range order is unlikely but correlation lengths should be large in very pure samples. 
The current patern in Fig. (6) can be observed by Bragg scattering of polarized 
neutrons or polarized x-rays. I estimate that with polarized neutrons the spin-flip 
cross-section in the circulating current phase at the (1, 1) Bragg peak is O(10~^) the 
nuclear cross-section. 

Another test would be evidence for local magnetic fields in regions 4 and 2 which 
are estimated to be 0(50 Gauss). As shown in Fig. (6), the local field is interstitial; 
there is no magnetic field either on Cu or O lattice sites. Muon spin resonance would 
be a way to look for interstitial fields if muons were to sit at the interstitial sites 
indicated in Fig. (6). 

As noted, the spectrum of current fiuctuations at q = as a function of a;, T and 
X is directly observable in Raman scattering. A direct test of the theory would be 
evidence for the fiuctuation spectra of Eq. (5.10) at large q and the difference in 
its projection on to the Cu and O sites obtainable by scattering in several Brillouin- 
zones mentioned above. In the section on NMR and inelastic neutron scattering, 
I mentioned that existing neutron scattering in Lai.85Sr.i5Cu04 is consistent with 
the magnetic fluctuations^^ derived from (5.10). More detailed tests, especially with 
scattering over a large range in momentum and frequency in YBa2Cu306.9 which 
shows no nesting related peaks in q-space, are suggested. 

In very pure samples of YBa2Cu306.7, and (248) which at stoichiometry behaves as 
an underdoped material, the sliver region 4 between regions 1 and 2 in Fig. (1) where 
the resistivity falls below the extrapolation from high temperatures clearly appears. 
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The cross-over to the circulating current phase in such samples should not be two 
broad. One should look for signatures of this in thermodynamic experiments, specific 
heat and magnetic susceptibility which should show a bump near the cross-over and 
a significant decrease below. Such samples at low temperatures would be particularly 
suitable to look for direct evidence discussed above the circulating current phase. 

The conjecture about magneto-transport^^ and the behavior of the critical fluc- 
tuations in a magnetic field can be tested by a Raman scattering experiment in a 
magnetic field. The polarized part of the spectra proportional to the magnetic field 
should acquire singular low energy, low temperature form. 

Some of the other tests of the theory have already been mentioned. These in- 
clude: (i) Improved angle resolved photoemission experiments to verify Eq. (6.4) 
and its modifications due to impurities deduced in Ref. (63). (ii) Measurement of 
the electronic specific heat in low Tc copper-oxides (for example the single layer Bi 
compound with Tc ^ lOK near ideal composition) to see the T£nT contribution to 
the electronic specific heat, (iii) Controlled single-particle tunneling experiments to 
see the spectrum of the "glue" for superconductivity. It should be mentioned that 
optical conductivity and Raman scattering experiments for a; > 2 A can also be used 
to deduce the spectrum of the "glue" . 
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Appendix A: One Electron Bandstructure and Eigen- 
vectors 

The bandstructure and the eigen- vectors for Hq of Eq. (2.2) for tpp/tpd << 1 are given 
here. For tpd = 0, the bonding "b" and the anti-bonding "a" bands have dispersions. 

6°b(k) ^ ± ((Ao/2)^ + Ul,sl^{k)f" (A.l) 

where s^y(k) = sin^ (^) the non-bonding band, "n" is non-dispersive with 
energy —Aq/2. The eigenvectors may be specified by the band annihilation operators 
in terms of these in the orbitals dj^^, Pxko-) Pyko- ^ ™ ^'i- (2-5) ■ 

The coefficients in (2.5a-b) for tpp = (specified by a superscript o) are 

<d(k) = (f + e°(k)) /Na(k) 

uL,y(k) =2itpd/N,(k) 

u°bd(k) =(t + eS(k))/Nb(k) 

u°bx,y(k) = -2itpd/Nb(k) 

where Na,b(k) = [(f e° b(k))^ + 4t^jS^ y(k)] ^2. The non-bonding orbital has an 
energy —A/ 2 and is annihilated by 

i(sin(kya/2)pxk - sin(kxa/2)pyk)/sxy(k) (A.3) 

The changes in the coefficients in (2.5a-b) are calculated to first order in tpp using 
as a perturbation the — hopping Hamiltonian 

H, = 4tpp E sin {^^^ sin (^^ j p+ ,Py\, + h.c. (A.4) 
u,d(k) =u°d(k)+fab(k)uSd(k) 

Uax(k) =s,(u°,(k)-f,b<y(k))+f,XySy (A.5) 

Uay(k) = Sy(u°y(k) - fabU^(k)) - facU°ySx 
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For the coefficients in (2.5b) replace a b in (A. 4). In (A. 4) 



fba(k) , . 

f n N 4tpp I 2 J I 2 J . kxa • kva 

fac(k) =^^rz^ .s?...fki sm- 



= -fbc(k) 
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Appendix B: Exchange Hamiltonian in r -space 

This appendix is primarily the work of Q. Si. Here the exchange Hamiltonian (3.18) 
is derived from the strong-coupling limit where the high energy states described in 
Sec. (3a) are eliminated by a canonical transformation. 

As usual in such a procedure, we write the Hamiltonian as 

H = Hiow + Hhigh + Hmix (B-1) 

where Hio^ contains the low energy states we wish to keep, i.e. states d]'";^|0 >, 
da^lO >■ Hhigh the states we wish to discard and Hmix connects states in Hio^ and 
Hhigh ■ We introduce a canonical transformation 

H = e^^H e~'^ (B.2) 

such that to linear order in Hmix, matrix elements connecting states of Hiow and Hhigh 
vanish. This requires that S be determined by 

+ i [S, Hio^ + Hhigh] = 0. (B.3) 

The transformed Hamiltonian, to second order in Hmix is 

H = Ho + i[SHmix], (B.4) 

is the kinetic energy (3.4) Hmix, which besides the allowed part in the low energy 
space (3.16) mixes states of Hio^ and Hhigh- We give the results for elements of M 
in the approximation that all the "Neglected stated' specified in Sec. (3a) except the 
zero-hole state (with energy 0) are assumed to be infinitely high compared to the 
low-energy states: the two one-hole states at energy =i=A — n and the two-hole state 
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at energy E,^ = V — 2/x: (No essential difference arises in tlie more general and messy 
situation.) Accordingly, we write the kinetic energy in terms of operators dii^ and 
d2iCT using (3.6), (3.7) and the first term of (3.17). 



div = ^sgn(70+d2i-a + di^0oi 
^ta = 75Sgn(T0,^dii_<, + d^i^0i. 



(B.5) 



The bare kinetic energy or Hmix, Eq. (3.4) is 



H • — V- • t^°d:^D- +t^?'^d:^"d- +tP°D:^D- 
+ h.c. 



(B.6) 



Inserting (A. 5) in (A. 6), we solve for S in Eq. (A. 3) by taking matrix elements 
between the states of Hiow + Hhigh of known energy. S is then inserted in (A. 4). The 
second term gives the exchange Hamiltonian of the form (3.19). The part in spin cr 
space is isotropic because of rotational invariance in <t space. The part in (di — d2) 
space is specified below 



y X z/ 



M has the form 



/ M 


V 



yy 







1 

y 

\< 




(B.7) 



(B.8) 



which we already discussed is the most general form allowable. 
We find that with E^ ± 2A = E± and E'^ + El^ = E'^ 



yy 
Ml 



= -2 t^D/E, + (tL) VE 

-XX = -2 t^d^DD/E^i - (t^D)^/E^ 

=(t:^d+t;^)/E.-(tL)VE 



(B.9) 



A rotation (3.32)-(3.33) about Ty which diagonalizes the kinetic energy to a-b 
space of Sec. (4), is used to get the Hamiltonian which on Fourier transforming gives 



74 



(4.2). For the special case A = 0, the rotation is by an angle j. In that case and if 
tdd — toD = 0, 

{J..-Jyy) =-tlJE^ (B.IO) 
iTzx — 2^dD/^</>- 

Then for (tz) = 1, A + 2B = 0. For the more general condition on A or 

X can always be found, so that the condition for a QCP derived in Sec. (4), i.e. 
±A + 2B + 4C = is fulfilled. 
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Figure Captions 



Fig. 1. Schematic Generic phase diagram of the quasi-two-dimensional Copper-Oxide 
compounds, x is the density of holes doped in the planes. Xc is the "optimum" 
composition. The antiferromagnetic phase and the superconducting phases, 
shown inside solid lines, occur through phase transitions. A series of cross-overs 
shown through dashed lines are discussed in the text. The impurity density, as 
inferred from the extrapolation of the high temperature resistivity to T = 0, 
i.e. assuming Mattheisen's rule decreases as x increases. The size of region 4 
decreases with increasing disorder for a given x. 

Fig. 2. The Zaanen-Sawatzky- Allen (ZSA) phase diagram for 3d-transition metal ox- 
ides, slightly modified and showing the schematic change in the position of the 
transition metal oxides going from left to right of the periodic table. The modifi- 
cation is that one of the axes is the ionic energy Ex defined through particle-hole 
spectra. ZSA used the charge transfer gap A defined through one-particle spec- 
tra. Ex < A due to particle-hole interactions. Ex/W is the more appropriate 
parameter to characterize the metal-insulator transition than A/W. 

Fig. 3. a Schematic one hole spectra (measured in photoemission) and one particle-spectra 
(measured in inverse photoemission) projected on to Cu states and oxygen states 
in Cu-0 compounds at 1/2-fiUing. 

Fig. 3.b Schematic one-hole spectra and one-particle spectra for a transition metal oxide 
far to the left of Cu say Cr-0, projected on to Cr and to O. 

Fig. 4. The bonding b and the antibonding band a for the two-dimensional band struc- 
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ture from one-electron theory in the hole representation with chemical poten- 
tial II. Under the q = transition discussed in the text, identical internal 
re-arrangements in each unit cell occur. So the band structure changes merely 
to the bands (3 and a shown. 

Fig. 5. (a) The calculated mean-field free-energy as a function of the real part of interband 
order parameter Tx described in the text. 

Fig. 5.(b) For a fixed Tx = T°, the free energy as a function of the imaginary part of the 
interband order parameter Ty. Ty takes the value for p > and a finite value 
below p < through a second order transition p is a parameter, defined in terms 
of the parameters of the Hamiltonian and for a given compound can be varied 
by varying the electron density, temperature or pressure. Ty 7^ corresponds 
to a circulating current pattern in the ground state shown in Fig. (10). 

Fig. 6. The deduced ground state current distribution pattern in the circulating current 
phase drawn for four cells. The -|- and — signs indicate magnetic fields pointing 
up and down. 

Fig. 7. The theoretical phase diagram in the pure limit. The effect of impurities is 
discussed in the text. 

Fig. 8. The unit cell of Cu-0 compounds in the x-y plane and the minimal orbital set: 
dx2_y2 orbital of Cu and a Px and a Py orbital of oxygen per unit cell. 

Fig. 9. (a) Exact representation of the interband susceptibility as a function of energy cu 
and momentum q. The lines are exact one-particle Green's functions and A is 
the complete (reducible) vertex. 
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(b) The interband susceptibility in the ladder diagram approximation to A of 

(a) . 

(c) Elementary self-energy and vertex corrections neglected in (b) . 

Fig. 10. (a) Interband absorption spectrum near the excitonic threshold in the approxima- 
tion that band b is infinitely massive - after Ref. (52). 

(b) Interband absorption spectrum near the excitonic ledge with finite hole 
mass. The excitonic ledge shifts to lower energy as the interband interactions 
increase. 

Fig. 11. Processes for analysis of the low energy Hamiltonian; (a), (b), (c) are processes 
considered for the Boson self-energy, (d) for the Fermion self-energy, (e) is the 
lowest order correction to the Boson-Fermion vertex, (f) Leading self-energy 
due to anharmonic interaction between the fluctuations. 

Fig. 12. Elementary processes for optical conductivity. 

Fig. 13. (a) The calculated optical conductivity based on the theory which reproduces re- 
sults based on phenomenology in Ref. (12a). The parameters used are ~ 2eV 
A = 0.5 and Uc — 1200K. A soft cutoff is used. 



Fig. 13. (b) Experimental results for the optical conductibility in the basal plane, from Ref. 
(58). 
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